Unambiguous State Discrimination 
of two density matrices 
in Quantum Information Theory 



Den Naturwissenschaftlichen Fakultaten 
der Friedrich-Alexander-Universitat Erlangen-Niirnberg 

zur 

Erlangung des Doktorgrades 



vorgelegt von 
Philippe Raynal 
aus Lyon, Frankreich 



ii 

Quantum Information Theory Group 
Theoretische Physik I 



Lehrstuhl fur Optik 
Institut fur Optik, Information und Photonik 
Max Planck Forschungsgruppe 



Erlangen 2006 



iii 



Als Dissertation genehmigt von den naturwissenschaftlichen Fakultaten der Universitat 
Erlangen-Niirnberg 



Tag der miindlichen Priifung : 16.08. 2006 

Vorsitzender der Promotionskommission: Prof. Dr. D. P. Hader 

Erstberichterstatter: Prof. Dr. N. Liitkenhaus 

Zweitberichterstatter: Prof. Dr. D. BruB 



iv 



Abstract 



Quantum state discrimination is a fundamental task in quantum information theory. The signals 
are usually nonorthogonal quantum states, which implies that they can not be perfectly distin- 
guished. One possible discrimination strategy is the so-called Unambiguous State Discrimination 
(USD) where the states are successfully identified with non-unit probability, but without error. 
The optimal USD measurement has been extensively studied in the case of pure states, especially 
for any pair of pure states. Recently, the problem of unambiguously discriminating mixed quan- 
tum states has attracted much attention. In the case of a pair of generic mixed states, no complete 
solution is known. In this thesis, we first present reduction theorems for optimal unambiguous 
discrimination of two generic density matrices. We show that this problem can be reduced to 
that of two density matrices that have the same rank r in a 2r-dimensional Hilbert space. These 
reduction theorems also allow us to reduce USD problems to simpler ones for which the solu- 
tion might be known. As an application, we consider the unambiguous comparison of n linearly 
independent pure states with a simple symmetry. Moreover, lower bounds on the optimal failure 
probability have been derived. For two mixed states they are given in terms of the fidelity. Here 
we give tighter bounds as well as necessary and sufficient conditions for two mixed states to 
reach these bounds. We also construct the corresponding optimal measurement. With this re- 
sult, we provide analytical solutions for unambiguously discriminating a class of generic mixed 
states. This goes beyond known results which are all reducible to some pure state case. We how- 
ever show that examples exist where the bounds cannot be reached. Next, we derive properties 
on the rank and the spectrum of an optimal USD measurement. This finally leads to a second 
class of exact solutions. Indeed we present the optimal failure probability as well as the optimal 
measurement for unambiguously discriminating any pair of geometrically uniform mixed states 
in four dimensions. This class of problems includes for example the discrimination of both the 
basis and the bit value mixed states in the BB84 QKD protocol with coherent states. 
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Zusammenfassung 



Quantenzustandsunterscheidung ist eine fundamentale Aufgabe der Quanteninformationstheo- 
rie. Die Signale sind normalerweise nicht-orthogonale Quantenzustande, d.h. sie konnen nicht 
perfekt unterschieden werden. Eine der moglichen Unterscheidungsstrategien ist die so genan- 
nte Eindeutige Zustandsunterschiedung (Unambiguous State Discrimination - USD), bei der die 
Zustande mit einer Wahrscheinlichkeit kleiner als eins erfolgreich erkannt werden, allerdings 
fehlerfrei. Optimale USD-Messungen fur reine Zustande sind ausfuhrlich untersucht worden, 
insbesondere fur jedes Paar von reinen Zustanden. Vor kurzem hat die Aufgabenstellung der 
eindeutigen Zustandsunterscheidung gemischter Zustande viel Aufmerksamkeit auf sich gezo- 
gen. Im Falle eines Paares von allgemeinen gemischten Zustanden ist keine vollstandige Losung 
bekannt. In dieser Doktorarbeit legen wir zuerst Reduktionstheoreme fur optimale eindeutige 
Unterscheidung von zwei allgemeinen Dichtematrizen vor. Wir zeigen, dass diese Aufgaben- 
stellung reduziert werden kann auf diejenige von zwei Matrizen, die denselben Rang r in einem 
2r-dimensionalen Hibert-Raum haben. Diese Reduktionstheoreme ermoglichen uns ebenfalls, 
USD-Aufgaben auf einfachere zuruckzufiihren, fur die die Losung moglicherweise bekannt ist. 
Der eindeutige Vergleich von n linear abhangigen reinen Zustanden mit einfacher Symmetric 
wird als Anwendung behandelt. Dariiber hinaus wurden untere Grenzen fur die optimale Fehler- 
wahrscheinlichkeit entwickelt. Fur zwei gemischte Zustande werden diese in Form der Fidelity 
angegeben. Hier geben wir engere Grenzen an, ebenso wie notwendige und ausreichende Be- 
dingungen fur zwei gemischte Zustande, diese Grenzen zu erreichen. Wir konstruieren ebenfalls 
die entsprechende optimale Messung. Zusammen mit diesem Ergebnis prasentieren wir ana- 
lytische Losungen fur die eindeutige Unterscheidung einer Kategorie allgemeiner gemischter 
Zustande. Dies geht iiber bekannte Ergebnisse hinaus, die alle auf reine Zustande zuriickfiihrbar 
sind. Wir zeigen allerdings, dass es Beispiele gibt, bei denen die Grenzen nicht erreicht werden 
konnen. Als nachstes leiten wir Eigenschaften des Rangs und des Spektrums einer optimalen 
USD-Messung her. Dies fiihrt schlieBlich zu einer zweiten Kategorie exakter Losungen. Wir 
zeigen die optimale Fehlerwahrscheinlichkeit auf, ebenso wie die optimale Messung, um jedes 
Paar geometrisch gleichformiger gemischter Zustande in vier Dimensionen zu unterscheiden. 
Diese Kategorie von Aufgabenstellungen schlieBt zum Beispiel die Unterscheidung von sowohl 
der basis- als auch der bit value-gemischten Zustande des BB84-QKD-Protokolls mit koharenten 
Zustanden ein. 
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Chapter 1 
Prologue 



Physics attempts to describe the world with the language of mathematics. Given a system an 
observer summarizes his knowledge in an abstract mathematical object, the so-called 'state'. At 
a given point in time this observer may decide to acquire information about the system. Such an 
acquisition of information is called a measurement. In that sense, Quantum Mechanics is con- 
cerned with knowledge, and the two pillars of Quantum Mechanics are states and measurements. 

Information Theory started in the late 1940's boosted by the second world war and its needs 
for communication and computational power. Information Theory addresses the fundamental 
questions of the transmission, processing and coding of information. 

It is therefore quite natural that Quantum Mechanics and Information Theory finally merge 
to describe the production, the transmission and the detection of information as well as its pro- 
cessing and coding. Quantum Information Theory was born. 

1 .1 Quantum Information Theory 

Since no information-theoretic formulatiorQis yet available, Quantum Information Theory (QIT) 
is formulated on the basis of four postulates that mathematically describe a physical system, its 
evolution and measurements that can be performed on it. Let us now review these four postulates 
EL 

Postulate 1 Hilbert space 

Associated to any isolated quantum system is a Hilbert space known as the state space of 
the system. The system is completely described by a unit vector | Y) called the state vector 
in the state space. 



See the work of R. Clifton, J. Bub and H. Halvorson or the work of A. Grinbaum for two appealing attempts. 
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Postulate 2 Unitary evolution 

The evolution of a closed (i.e. an isolated system having no interaction with the environment) 
quantum system is described by a unitary transformation. That is, if | Y) is the state at time 
t, and |Y ; ) is the state at time t', then |Y ; ) = E/|Y) for some unitary operator U which 
depends only on t and t' . 

Postulate 3 Measurement 

A measurement is described by a collection {M m } of measurement operators. These op- 
erators are acting on the state space of the system being measured. The index m refers to 
the measurement outcomes that may occur in the experiment. If the state of the quantum 
system is | Y) immediately before the measurement then the probability that result m occurs 
is given by p(m) = (Y|M^M m |Y), and the state of the system after the measurement is 

M„,|Y) 
y/C¥\MlM m \T)' 

Moreover the measurement operators satisfy the completeness equation, L m M^Mm = 1. 

Note that in Quantum Information Theory the measurement operators {M m } are often called 
Kraus operators 0. 

Postulate 4 Composite system 

The state space of a composite quantum system is the tensor product of the state spaces 
of the component quantum systems. That is, if we have systems numbered 1 through n, 
and system number i is prepared in the state |Y,), then the joint state of the total system is 
|Yi)®|Y 2 >®"-®|Y B >. 

Note that, unlike in Quantum Mechanics, observables do not have a crucial role in Quantum 
Information Theory. Moreover, in general, we can consider the state of a system to be not only a 
vector state but a classical mixture of vector states. The notion of density matrices then is useful 
as we will see in the next subsection. Measurements are the core of Quantum Information theory 
because it is through a measurement that we learn information about a system. Therefore, we 
also introduce the mathematical language used to described a measurement. 

1.1.1 Ensemble of quantum states and density matrix 

Let us suppose a quantum system is in the state |Y,-) chosen in a set of states {|Y,)}. We can 
imagine that the appearance probabilities r/,- of each state of the set are in general different. 
We then summarize our knowledge on the system with the ensemble {|Y,-), 77,-}. It is called an 
ensemble of the system. If the ensemble is composed of only one state (and of course its a priori 
probability equals 1), the state is called pure. If not, one speaks of mixed states that is to say a 
classical mixture of pure states. To efficiently describe a mixed state, we use an operator instead 
of a vector state, the so-called density matrix. 



1.1. Quantum Information Theory 
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Definition 1 Density matrix 

Let us consider a system with ensemble {|Y;),TJ;}. The state of the system can then be described 
in a compact form by the density matrix 

p = En/|Y,-)<Y,-|. (i.i) 

i 

Such a density matrix possesses the three important properties 

Tr(p) = 1 (Normalization), (1.2) 
p > (Positivity), (1.3) 
Tr(p 2 ) = 1 : p = |Y)(Y| (Purity), (1.4) 

where > means positive semi-definite. Actually, the state ensemble of a system is not unique. 

Theorem 1 Unitary freedom in the state ensemble 

The sets {|Y;),T},-} and V,-} generate the same density matrix if and only if there exists a 

unitary transformation U such that 

VVi^^Uij^njiYi). (i.5) 

j 

Equivalently, 

Corollary 1 Unitary freedom in the state ensemble of a density matrix 

The two density matrices L/ty |%')(%'| and Y,i v i\*&i)( < &i\ describe the same state if and only if 
there exists a unitary transformation U such that 

y/Vt\<S>i)=ZUij^rj\Yi). (1.6) 

j 

1.1.2 Generalized measurements - POVM 

The third postulate of QIT, and its measurement operators E m , can be used to define the posi- 
tive semi-definite operators E m = M\ x M m . The set {E m } m is called a Positive Operator- Valued 
Measure (POVM) El HI and each operator E m , a POVM element. On one hand, the fact 
that the probabilities p(m) = (Y|.E m |Y) are real and positive is expressed by the positivity of 
the POVM elements {E m } m . On the other hand, the fact that probabilities add up to one is ex- 
pressed by the completeness relation Y, m E m — 1- Indeed, the sum of the probability p(m) is 
^ m p(m) = L m (Y|£ m |Y) = (Y|I OT £ m |Y) = (Y|Y) = 1. An important property of a POVM 
element is that its spectrum is upper bounded by 1 . Otherwise, it is clear that the expectation 
value (Y\E m \W) would exceed unity which contradicts the requirement that a probability is less 
than 1 . We finally give a general definition of a POVM. 
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Definition 2 POVM 

A Positive Operator-Valued Measurement (POVM) is a set of positive semi-definite operators 
{E m }m such that 

E k > (Positivity) (1.7) 
= 1 (Completeness relation) (1.8) 

k 

The probability to obtain the outcome kfor a given state pi is then given by 

p{k\i) = Tr{E kPi ). (1.9) 

In the previous formula, Tr(.) stands for the trace. A POVM is also called a generalized 
measurement since it is the most general description of a measurement. Indeed, projective mea- 
surements, usually encountered in Quantum Mechanics are, in the above formalism, merely a 
special case where E m E n = 8 mn , E^ = E m . Such a projective measurement is called a Projection 
Valued Measure (PVM). Nevertheless, a generalized measurement can also be described by a 
projective measurement on an enlarged Hilbert space. A generalized measurement is then seen 
as a special case of projective measurements. The two pictures finally are equivalent as long as 
the Hilbert space is not fixed. This is made precise in the following theorem due to Naimark 

BH. 

Theorem 2 Naimark 's extension 

Given {E^} a POVM on a Hilbert space 3if, it exists an embedding of 3f into a larger Hilbert 
space 3$T such that the measure can be described by projections onto orthogonal subspaces in 
J?T. That is, there exist a Hilbert space 3¥ , an embedding <§ such that S3^ = 3f and a PVM 
{Rk} in 3t ', such that with P the projection defined by P3¥ = 3if, E^ = PRkP, Vfc. 

1.1.3 Definitions and notations 

Here we briefly fix some notations. Throughout this thesis, we will make an extensive use of the 
support S^p := support(P) of a Hermitian operator P. The support of a Hermitian operator is 
defined as the subspace spanned by its eigenvectors. We can moreover define the kernel 3tp := 
kernel (P) of a Hermitian operator P as the subspace orthogonal to its support. We also denote 
rp := rank(P) = dim {S^p), the rank of P. 

Next we define in a Hilbert space 3t? the sum and the intersection of two Hilbert subspaces 
3f\ and The sum 31C\ + 3ffi of the subspaces 31C\ and 3^2 is defined to be the set consisting 
of all sums of the form a\ + a%, where a\ E 3^\ and a2 G 3fii. 3f\ + 3?2 is a Hilbert subspace 
of 3f. The intersection 3^\ fl 3fii is defined to be the set consisting of all the elements a, where 
a E 3t\ and a E 3^i- H 3fii is a Hilbert subspace of 3^ . The complementary orthogonal 
subspace ( or orthogonal complement) of a subspace 5? in 3¥ ', written 5?^, is the set of all the 
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elements of orthogonal to 5? with respect to the usual euclidean inner product. We then have 
ffl = y © J? , the direct sum of the two orthogonal subspaces. Note that we use indifferently 
the notation J^p or S^p for a Hermitian operator P. 

We need to define a positive semi-definite operator. A Hermitian operator A acting on M' is 
positive semi-definite if and only if (Y|A|Y) > 0, for all |Y) in M 1 . In other words, a Hermitian 
operator is positive if and only if all its eigenvalues are positive or zero. We use the notation A > 
to say that an operator A is positive semi-definite. For such a positive semi-definite operator 
A. We can define its unique square root VA and decompose it into the form A = MM^ with 
M = VaU, for any unitary matrix U. Since the states p, and the POVM elements are positive 
semi-definite operators, we can introduce their square root and use the previous decomposition. 



1 .2 Unambiguous Quantum State Discrimination 

A quantum state describes what we know about a quantum system. Given a single copy of a 
quantum system which can be prepared in several known quantum states, our aim is to determine 
in which state the system is. This can be well understood in a communication context where 
only a single copy of the system is given and only a single shot-measurement is performed. This 
is in contrast with usual experiments in physics where many copies of a system are measured 
to get the probability distribution of the system. In quantum state discrimination (see for a 
review of quantum state discrimination), no statistics is built since only a single-shot measure- 
ment is performed on a single copy of the system. Actually there are fundamental limitations to 
the precision with which the state of the system can be determined with a single measurement. 
Whenever the possible quantum states are nonorthogonal, perfect discrimination of the states 
becomes impossible. This can be understood from the intuition that two non-orthogonal states 
have some probability to behave the same way. More precisely, if a quantum system is prepared 
in one of the two state |Y) and which are neither identical nor orthogonal, there is no mea- 
surement that perfectly determines in which state the system is. Mathematically, a measurement, 
that perfectly determines in which state the system is, is composed of two outcomes (i.e. two 
POVM elements) £V and that identify |Y) and |<I>) respectively with no errors. This means, 
in terms of probabilities, that 

(Y|£ Y |Y) = 1, (1.10) 

(3>|£*|<E>) = 1, (l.H) 
(Y|£ |Y) = 0, (1.12) 

(<J>|£y|<£> = 0. (1.13) 
If we express |<I>) in the basis {|Y), (Y^)}, Eqn. (1.1 1) becomes 

((Y|0)*(Y| + (Y ± |0)*(Y ± |)£ 4 ,((Y|<1))|Y) + (Y^O)^)) = 1 (1.14) 
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where * stands for complex conjugation. With the help of Eqn.(1.12) which is equivalent to 
E&\*¥) = since > (see proof in Appendix A), we obtain 

|(d)|Y- L )| 2 (Y- L |£ci>|Y- L ) = 1. (1.15) 

Since the spectrum of E& is upper bounded by 1, (Y- L |£ , 4 ) |Y- L ) < 1 and Eqn. dl.15b is fulfilled 
only if KOjY^)! 2 = 1 which contradicts the assumption that |Y) and |<I>) are non-orthogonal. 

The immediate consequence of this limited precision is to resort to various state discrimina- 
tion strategies depending on what one really wants to learn about the state. Given a strategy, we 
finally have to optimize the measurement with respect to some criteria. 



A/ice 




Bob 


Po 
Pi 


signaL 


Measurement 



Figure 1.1: Two parties Alice and Bob want to communicate 



The basic scenario involves two parties Alive and Bob who want to communicate (see 
Fig. Il.lb . Alice prepares a quantum system in a state, member of a set of states known by 
Bob. In general Alice does not prepare each state with the same probability. We speak of an a 
priori probability. She sends a quantum system to Bob who performs a measurement in order to 
obtain the information he wants. In other words, a state ensemble of a quantum system is given 
and we want to determine the state of that system. In his famous book published in 1976 J3j, 
Helstrom established the mathematical bases of such detection tasks. He introduced the notion of 
Bayes' cost function which can describe any discrimination strategy. The idea is the following. 
For each possible outcome k conditioned on a signal state j, a price to pay Q/ is associated. If 
Ckj is positive, Bob has to pay Alice. If Q 7 is negative, Bob earns money. To set up a strategy 
corresponds to give the Bayes' cost matrix Ckj. Related to this matrix, the Bayes' cost function, 
given by 

C = Y,ljCkjp(k\j), (1.16) 

kj 

represents the total price that Bob has to pay to Alice. Information about a state is represented 
by an outcome k conditioned on a signal state j. It then appears clear that, depending on which 
information really matters to Bob and Alice, the strategy or, equivalently, the Bayes ' cost matrix 
Ckj will change. The aim for Bob is of course to minimize the prize he has to pay to Alice. To 
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minimize the Bayes' cost function C while the a priori probability r\j and the states p ; - are fixed, 
Bob is only free to change his measurement. In this thesis, we play the role of Bob who wants to 
find the optimal measurement to lose a minimal amount of money. 

The Bayes ' cost matrix Ckj depends on the strategy adopted by Alice and Bob. For instance, 
Bob might want to know which state was sent with the minimum error probability. This strategy 
is called Minimum Error Discrimination (MED) [3] - see Fig. 1 1 .2.1 In MED, the Bayes' cost 
matrix Ckj is given by 

<>■ {Hi. 



A/ice 




Bob 


Po 


s/gnaL 


Po 


Pi 


Pi 



Figure 1.2: Two possible outcomes in the scenario of Minimum Error Discrimination 



Alternatively, one might consider an error-free discrimination of the signal states. In this 
strategy, the measurement can either correctly identify the state or send out a flag stating that 
it failed to identify the state. A correct identification of the state is called a conclusive result 
while a failure to identify the state is known as an inconclusive result usually denoted by '?' 
or 'don't know'. The objective then is to minimize the probability of inconclusive result, the 
so-called failure probability. This strategy is called Unambiguous State Discrimination (USD) - 
see Fig. II .31 The coefficients of the non square (j = 0, 1 and k = 0, 1,?) are Bayes' cost matrix 
C k j are 

r o k=j, 

C kj =l 1 k=?,Vj, (1.18) 
^ oo otherwise. 

Note that the coefficients ; where k, j = 0, 1 are set to infinity in order to impose the error- free 
conditions p(k\j ^ k) = 0, k,j = 0, 1 to obtain a non diverging Bayes' cost function. 

We can list another task related to state discrimination where we are given a finite number of 
identical copies of an unknown state in a d-dimensional Hilbert space. Our goal is to estimate 
the actual state with the maximum accuracy, which is often quantified by the fidelity between 
the actual state and the estimated state (see chapter 2 for a definition of the fidelity). Since the 



8 



1. Prologue 



Alice 



Bob 



Po 
Pi 



Po 

Pi 
? 



Figure 1.3: Three possible outcomes in the scenario of Unambiguous State Discrimination 

state to estimate can be any state in the J-dimensional Hilbert space, one has to average the 
accuracy over all the possible states of the J-dimensional Hilbert space. This scenario is known 
as Quantum State Estimation JEHU (see Ref. OH [TQ CE2 for other scenarios). 

Let us add another comment. The fact that non-orthogonal quantum states are not perfectly 
distinguishable also has benefits. It leads in particular to secure Quantum Key Distribution 
(QKD) in a cryptographic context |fT3ll . The security in classical computer science is ensure by 
the complexity of some task like factorization of big prime numbers. In QKD, the security is due 
to the quantum laws of Nature and does not anymore rely on the assumption of eavesdropper's 
limited computational power. 

In general, the optimal measurements for a given strategy depends on the quantum states and 
the a priori probability of their appearance. For a given strategy and a given state ensemble, the 
task is to find the measurement which minimizes the Bayes' cost function. Such a measurement 
(it might not be unique) is called an optimal measurement. 

In this thesis, we are interested in the unambiguous discrimination of two known mixed 
quantum states. Therefore the task is to find an optimal measurement that minimizes the failure 
probability. The problem of unambiguously discriminating pure states with equal a priori 
probabilities was formulated in 1987 by Dieks 031 and Ivanovic lfT5l and elegantly solved 
by Peres ifToll . Seven years later, Jaeger and Shimony presented the general solution for two 
pure states with different a priori probabilities ifTTl . Shortly after this result, Chefles and 
Barnett showed that only linearly independent pure states can be unambiguously discriminated 
[fT8l . Finally Chefles provided the optimal failure probability and its corresponding optimal 
measurement in the case of n symmetric states [19]. The enumeration of analytical results 
for USD of pure states scenarios already ends here even if an algorithm for the case of three 
pure states was proposed by Peres and Terno in 1998 ll20l . In fact, since Sun's work in 2002 
[|2Tll22l . it is known that USD (of both pure and mixed states) is a convex optimization problem 
Il23ll24ll25l . Mathematically, this means that the quantity to optimize as well as the constraints 
on the unknowns are convex functions. Practically, this means that the optimal solution can be 
extremely efficiently computed. This is therefore a very useful tool. Nevertheless our aim is to 
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understand the structure of USD, to relate it to neat and relevant quantities and to find analytical 
solutions. 

The case of mixed states recently attracted more attention. But until this present work, no 
optimal measurements for mixed states has been found unless the USD problem can be reduced 
to some known pure state case. This reduction comes from simple geometrical considerations 
and can be summarized in three theorems. Important examples of such reducible problems 
are Unambiguous State Discrimination of two mixed states with one -dimensional kernel [|26l . 
Unambiguous State Comparison [|27l |28l |29l (see Ref. Il27l l30l I3T1 for the unambiguous 
comparison of unknown states), State Filtering (32l [33l [34l and Unambiguous Discrimination of 
two subspaces ll35l . This four cases are all reducible to some pure state case and can therefore 
be solved. To specify that a USD problem is not reducible by means of our three reduction 
theorems, we use the expression 'USD of generic density matrices'. Lower and upper bounds on 
the failure probability to unambiguously discriminate two density matrices are also known. In 
2004, Eldar derived necessary and sufficient conditions for the optimality of a USD POVM |[36l . 
Unfortunately these conditions appear rather difficult to solve. In contrast to the MED problem, 
which is already solved for any pair of mixed states AH [371, the optimal USD of mixed states is 
an open problem. 



1.3 Results 

We outline here the six main results derived in this thesis. 

1) Three reduction theorems to reduce the dimension of a USD problem 

2) Unambiguous comparison of n pure states with a simple symmetry 

3) First class of exact solutions 

4) Second class of exact solutions 

5) A fourth, incomplete, reduction theorem 

6) USD and BB84-type QKD protocol 
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Three reduction theorems to reduce the dimension of a USD problem [Chapter 3] 

As seen in the previous section, only few analytical optimal solutions in Unambiguous State 
Discrimination are known. For pure states scenarios, only two classes of exact solutions have 
been provided so far. They are the solutions for USD of two pure states ifTTl and USD of n 
linearly independent symmetric pure states lfT9ll . In the case of mixed states, there are actually 
four known solutions: unambiguous discrimination of two mixed states with one -dimensional 
kernel ll26l . unambiguous comparison of two pure states Il27ll28ll29l . state filtering Il32ll33l l34ll 
and unambiguous discrimination of two subspaces ll35l . It seems surprising that research on 
USD of pure states has been less successful than work on USD of mixed states! A solution to 
this apparent paradox is given by our first result. Indeed these four optimal solutions in USD of 
mixed states only require the optimal solution for USD of two pure states. More generally, we 
prove that the problem of discriminating any two density matrices can be reduced to the problem 
of discriminating two density matrices of the same rank r in a 2r-dimensional Hilbert space. 
This introduces the notion of standard USD problem. Such a standard USD problem is proposed 
as a starting point for any further theoretical investigation on USD. That way, we can avoid to 
deal with trivial or already known classes of solutions. The reductions are of three types and 
can be summarized in three theorems. In few words, the reduction theorems work as follows. 
In a first reduction theorem, we split off any common subspace between the supports of the two 
density matrices po and p\. In a second reduction theorem, we eliminate, if present, the part of 
the support of pi which is orthogonal to the support of po and vice versa. In a third reduction 
theorem, if two density matrices are block diagonal, we decompose the global USD problem into 
decoupled unambiguous discrimination tasks on each block. 

Unambiguous comparison of n pure states with a simple symmetry [Chapter 3] 

We are given n pure quantum states {|Y,)} which occur with a priori probabilities {pi}. We 
would like to know without error whether these states are all identical or not. Actually the task 
of unambiguously comparing any two pure states can be elegantly solved by use of the second 
and third reduction theorems, as Kleinmann et al. showed in il28ll . Stimulated by their idea, we 
investigate the case of n pure states having some simple symmetry. In fact we prove that the 
comparison of n linearly independent pure states with equal a priori probabilities and equal and 
real overlaps can be reduced to n unambiguous discriminations of two pure states and then be 
solved. The question to know whether any unambiguous comparison of pure states is always 
reducible to some pure state cases remains opened. Let us add here that, as Kleinmann et al. 
indicated in ll28l . the unambiguous comparison of mixed states is generally not reducible to 
some pure states case. 

In this thesis, we provide two classes of exact solutions for unambiguously discriminating 
two generic density matrices. These two classes are the only two classes known until now. 



1.3. Results 
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First class of exact solutions [Chapter 4] 

We consider the problem of unambiguously discriminating two density matrices po and pi with a 
priori probabilities r\o and r\\. We define the fidelity of the two states as F = Tr( */ y/popi y/po) ■ 
We provide three lower bounds on the failure probability in three regimes of the ratio between 



< 



the a priori probabilities defined as ^ < Tr( ^ p0 \ 1^^£°± < JlK < and 



Tr(PoPl) 

y^- For each regime, we give necessary and sufficient conditions for the failure probability of 
unambiguously discriminating two mixed states to reach the bound. With that result, we give the 
optimal USD POVM of a wide class of pairs of mixed states. This class corresponds to pairs of 
mixed states for which the lower bound on the failure probability is saturated. This is the first 
analytical solution for unambiguous discrimination of generic mixed states. This goes beyond 
known results which are all reducible to some pure state case. Note that any pair of mixed state 
does not saturate the bounds. The necessary and sufficient conditions take the simple form of 
the positivity of the two operators po — Ct \/y/poPiy/po and Pi — ^ \/ y/PiPoy/Pi where a equals 
Tr (APo) ^ fq± an( j _^F_^ m g rs ^ seconc i anc [ third regime, respectively. 

Second class of exact solutions [Chapter 5] 

We derive a second class of exact solutions. This class corresponds to any pair of geometrically 
uniform mixed states without overlapping supports in a four dimensional Hilbert space. In short, 
two geometrically uniform mixed states are two unitary similar density matrices po and pi = 
U PqU where the unitary matrix U is an involution i.e. U 2 = 1. We find that only three options 
for the expression of the failure probability exist. First, if the operators po — \J^Fq and p\ — 

J^F\ are positive semi-definite, then the pair of density matrices falls in the first class of exact 

solutions. If this is not the case, either the operator PqUPq has one positive and one negative 
eigenvalue or it has two eigenvalues of the same sign. In the former case, we can give the optimal 
failure probability in terms of the eigenvalues and eigenvectors of PqUPq. In the later case, no 
unambiguous discrimination is possible and the failure probability simply equals unity. For these 
three cases, we provide the optimal failure probability as well as the optimal measurement. 

A fourth, incomplete, reduction theorem [Chapter 5] 

The two USD POVM elements Eq and E\ have a rank less or equal to the rank of and 
respectively. This defines the notion of maximum rank of Eq and E\. We establish a 

theorem stating that if the two operators po — <\J^Fq an d Pi — \J~^F\ are not positive semi- 
definite then the two USD POVM elements Eq and E\ can not have both maximum rank. A 
corollary can be derived assuming a standard USD problem. In that case, if the two operators 
Po — y^^b an d Pi — y^i are not positive semi-definite then there exist one eigenvector of 



12 1. Prologue 

E-7 with eigenvalue 1 and one eigenvector of either Eq or E\ with eigenvalue 1, too. From the 
completeness relation fulfilled by the measurement operators, it follows that we can split off the 
two-dimensional subspace spanned by these two eigenvectors from the original USD problem. 
This could lead to a fourth reduction theorem. 'Could' because it remains to fully characterize 
these two eigenvectors cited above. So far, we can only prove their existence. If one could 
characterize them, a way to solve analytically any USD problem would be available. Indeed, we 
start from a general USD problem of two mixed states. We use the three first reduction theorems 
to bring it to standard form. We then check the positivity of the two operators po — \J~^Fq an d 

Pi — yF^F\. If the positivity is confirmed, then the pair of density matrices falls in the first 
class of exact solutions. If the two operators are not positive, we can use the fourth reduction 
theorem to get rid of two dimensions corresponding to the two eigenvectors mentioned above. 



At that point, we check the positivity of the two operators Pn — \ ^f^n and p[ — \ ^j-F! of the 

v "o V "i 



reduced problem. We see here a constructive way to solve any USD problem. If the two operators 

Pq — J^tFq and p[ — y^j^i °f me reduced problems never turn out to be positive, we end up 

with only two pure states and we can therefore always find the optimal measurement. The full 
characterization of the two eigenvectors involved in this incomplete reduction theorem is of great 
importance. 



USD and BB84-type QKD protocol [Chapter 6] 

The Bennett and Brassard 1984 cryptographic protocol ll38l provides a method to distribute 
a private key between two parties and allow an unconditionally secure communication. We 
consider in this thesis the implementation of a BB84-type QKD protocol that uses weak 
coherent pulses with a phase reference ll39l . In that context, two important questions related 
to unambiguous state discrimination can be addressed. First, 'With what probability can an 
eavesdropper unambiguously distinguish the basis of the signal?' and second 'With what 
probability can an eavesdropper unambiguously determine which bit value is sent without 
being interested in the knowledge of the basis?' These two questions can be translated in 
some unambiguous discrimination task of two geometrically uniform mixed states in a four 
dimensional Hilbert space. We answer these two questions providing useful insights for further 
investigations on practical implementations of Quantum Key Distribution protocols. 

The structure of this thesis is the following. In chapter 2, we mathematically define the prob- 
lem of USD. We then review the known results on unambiguous discrimination: unambiguous 
discrimination two pure states, unambiguous discrimination of n symmetric states and a few gen- 
eral properties. In chapter 3, we present our three reduction theorems. They allow us to solve 
special tasks in quantum information theory such as, e.g. state filtering, unambiguous discrimi- 
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nation of two pure states, unambiguous discrimination of n pure states with a simple symmetry 
and unambiguous discrimination of two subspaces. All these tasks are related to the unambigu- 
ous discrimination of two mixed states which can be reduced to the unambiguous discrimination 
of some pure states only. We also define a standard form as a starting point for further inves- 
tigations in USD. In chapter 4, we derive lower bounds on the failure probability Q as well as 
necessary and sufficient conditions for the failure probability to reach those bounds. This pro- 
vides a first class of exact solutions for unambiguous discrimination of two generic mixed states. 
This class corresponds to pairs of mixed states for which the lower bound (one for each of the 
three regimes depending on the ratio between the a priori probabilities) on the failure probability 
Q is saturated. For this class we give the corresponding optimal USD measurement. In chapter 
5, we derive a fourth, incomplete, reduction theorem which, together with the first three reduc- 
tion theorems aims to solve in a constructive way any USD problem of two density matrices. 
Moreover we derive a second class of exact solutions. This class corresponds to any pair of two 
geometrically uniform states in four dimensions. In chapter 6, we give two examples of such 
an unambiguous discrimination of two geometrically uniform states in four dimensions. These 
examples are related to the implementation of the Bennett and Brassard 1984 cryptographic pro- 
tocol. In the last chapter, we summarize our results and propose directions for further research 
on USD of two density matrices. 
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Chapter 2 

Optimal Unambiguous State 
Discrimination 

The optimal USD measurement is known for two pure-state cases. On one hand, the optimal 
failure probability as well as the corresponding optimal measurement were provided by Jaeger 
and Shimony for any pair of two pure states with arbitrary a priori probabilities ifTTl . On the other 
hand, Chefles found the optimal failure probability and the corresponding optimal measurement 
for unambiguously discriminating n linearly independent symmetric pure states |fT9l . We present 
the basic properties of a USD measurement before reviewing the solution to these two pure-state 
scenarios. 



2.1 The USD measurement 

We consider a set of n G IN known quantum states {pi}, i = 1, with their a priori probabil- 
ities {rji}. We are looking for a measurement that either identifies a state uniquely (conclusive 
result) or fails to identify it (inconclusive result). The goal is to minimize the probability of 
inconclusive result. The measurements involved are typically generalized measurements de- 
scribed by a POVM which consists in a set of positive semi-definite operators {E{\ that satisfies 
the completeness relation J^E^ = 1 on the Hilbert space spanned by the states. The probability 
to obtain the outcome k for a given signal p ; - is then given by p(k\i) = Tr(pjE k ). We will of- 
ten refer to the states of the quantum system as signal states or even signals. This comes from 
the context of communication where the possible states of a quantum system correspond to the 
different signals sent to communicate. 

Let us now mathematically define what an Unambiguous State Discrimination Measurement 
is, its corresponding failure probability, and the notion of optimality. 
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Definition 3 A measurement described by a POVM {E k } is called an Unambiguous State Dis- 
crimination Measurement (USDM) on a set of states {pi} if and only if the following conditions 
are satisfied: 

• The POVM contains the elements {E? ,Ei,...E n } where n is the number of different signals 
in the set of states. The element £V is connected to an inconclusive result, while the other 
elements E lt i = 1, ..,n , correspond to an identification of the state p,. 

• No states are wrongly identified, that is Tr (pjE k ) =0 \/i^k i,k = 1, ...,n. 

Each USD Measurement gives rise to a failure probability, that is, the rate of inconclusive results. 
This can be calculated as 

Q[{E k }} :=£T7 ( Tr( P/ £ ? ). (2.1) 

i 

Definition 4 A measurement described by a POVM {E? p } is called an Optimal Unambiguous 
State Discrimination Measurement (OptUSDM) on a set of states {pi} with the corresponding a 
priori probabilities {fy} if and only if the following conditions are satisfied 

• The POVM {E° k pt } is a USD measurement on {p,} 

• The probability of inconclusive results is minimal, that is Q^E^}] = mmQ[{E k }\ where 
the minimum is taken over all USDM. 

Unambiguous state discrimination is an error-free discrimination. This implies a strong con- 
straint on the measurement. The fact that the outcome E k can only be triggered by the state p k 
implies that the support of E k is orthogonal to the supports of all the mixed states other than p k . 
This is a strong constraint for any USD measurement, not only the optimal one. To see that fact 
rigorously we need the following lemma. 

Lemma 1 For any positive semi-definite operators A andB, Tr (AB) = if and only if the support 
of the two positive semi-definite operators are orthogonal 

Tr(Afl) = S A _L S B . (2.2) 

Since a USD POVM satisfies Tx(E k pi) = Tr(E k p k )8ki to be an error-free measurement, a 
corollary of Lemma [Dean be derived. 

Corollary 2 A USD measurement described by the POVM {E k } on n density matrices {p,} is 
such that 

S Ek -L Sp ¥k ,\/i,k=l,...,n. (2.3) 

USD measurements are very sensitive in the sense that a small variation of a mixed state 
overthrows completely the error-free character of the already existing measurement. This is true 
for any USD measurement, not only the optimal ones. Let us now prove Lemma 1. 



2.2. Solution for two pure states 
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Proof of Lemma 1 If A and B are positive semi-definite operators, they are diagonalizable with 
eigenvalues oc,- > (i = \,...,rank(A)) and /3 7 > (j = \,...,rank{B)). Thus 

Tr(A5) = TrfXoftlT/XT/lE^I^X*,!) 

= £<*0y|<T f |4v)| 2 (2.4) 

'7 

vanishes if and only if {|4>,-)} and {|Y ; )} span orthogonal subspaces. ■ 



2.2 Solution for two pure states 



In the simple case of two pure states |Yo) and |Yi) with arbitrary a priori probabilities t\q and 
r/i, the optimal failure probabilities (see Fig. 12.11) to unambiguously discriminate them is given 
by 



< 2 opt = 7 ?1 + T ] o|(Yo|Y 1 )| 2 for J^< |(Y |Yi)|, 

V 

2 opt = 2 V ^77T|(Y |Y 1 )| for |(Yo|Yi)| < < 



e° pt = T?o + T?i|<Y |Yi)| 2 for 



rjo ~~ |(Y |Yi)|' 
V rio 



l(Y |Y 



(2.5) 
(2.6) 
(2.7) 



This result was derived by Jaeger and Shimony in 1995. When the two a priori probabilities are 
equal, it reduces to the well known equation 



g opt = |(Yo|Yi] 



(2.8) 



This solution is known as the Ivanovic-Diesk-Peres (IDP) limit since 1988. 

The optimal measurement (see Fig. 12.21) that realizes these optimal failure probabilities is 
given by 



Eb=|Yj-)(Yf- 
#1 = 

Ei = |Yi)(Yi| 




fbrW^<|(T |Ti; 



(2.9) 



Eq 

Ei 
E? 



KYo|Y, 



Yf)(Y 



IWIYo)| 2 
l-/f KYolYQI 
l(Y^|Ti>|2 



1 /\ M 



for |(Y |Yi)| < 




< 



|(Y |Yi 



(2.10) 
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E = 

Et = \Vi)(ri\ for 
£ ? = |Y )(Y | 



< , / 5i. 



l(Yo|Ti>| -V?]o" 



(2.11) 



i 



(*ol*i) 







Q opL ?7i I %|{*o|*i) 2 




(*o|*i 



■ m 



Q opt ='W + '//i|{*n|*i)f 



Figure 2.1: Optimal failure probability for USD of two pure states 




Figure 2.2: Basis vectors (Y] 1 ), |Yq ) and |?) of the three POVM elements £o> £i an d £? f° r 
the optimal USD measurement of two pure states when (Yq|Yi) > and |(Yo|Yi)| < < 



2.3. Solution for n symmetric pure states 
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2.3 Solution for n symmetric pure states 

Unambiguous discrimination can be consider for more than two states. The only requirement for 
an error-free discrimination is the linearly independence of the signal states as Chefles showed 
in 1998. An exact solutions can even be provided if the n quantum states happen to be symmet- 
ric. Symmetric states are states that can be written in terms of a generator |Yo) and a unitary 
transformation U such that U" = 1. The complete set of symmetric states can be written as 

|Y 7 ) = U\Wj-i) = U j \Y ), j=\,...,n-\ (2.12) 
|Y ) = tf|Y n _i), U n = l. (2.13) 

Note that we choose the a priori probabilities to be equal in order not to break the symmetry. For 
such symmetric states, we can introduced a suitable orthonormal basis {|%)}jt such that |Y ; ) = 
LlZlc k e 2in ^\yk) with Zk\c k \ 2 = 1 and U = l! l k z} Q e 2l <\y k ){y k \ ED. Note that the coefficients 

c k can be calculated thanks to the formula \c k \ 2 = -tL/,/ ■e 2wc *~ (Y y/|Y '_,•) . We define c ml -„ = 
min k c k and the optimal failure probabilities to unambiguously discriminate n symmetric states is 
then given by 

Q opt = n\c min \ 2 . (2.14) 

On the analytical side, some general properties of USD of mixed states were recently de- 
rived. We give here an overview of these results. First, there are the very general necessary and 
sufficient conditions for the optimality of a USD measurement derived by Eldar in 11361 . Unfor- 
tunately those conditions are pretty hard to solve. They can nevertheless be used to check the 
optimality of some USD POVM or, as we will do in chapter 5, to derive a new class of exact 
solutions. This class correspond to pairs of two Geometrically Uniform density matrices in four 
dimensions. Another general result on USD of two mixed states is the derivation of lower and 
upper bounds on the optimal failure probability. The lower bounds are expressed in terms of the 
fidelity. Therefore we first introduce this quantity. The upper bound is presented in term of the 
failure probabilities of some pure state case. 

2.4 Necessary and sufficient conditions for the optimal- 
ity of a USD measurement 

Necessary and sufficient conditions for an optimal measurement that minimizes the probability 
of inconclusive result can be derived using argument of duality in vector space optimization ll36l . 
These conditions are valid for any number of mixed states. Let us now state the theorem. 

Theorem 3 Let {pi}, I <i<n denote a set of density operators with their a priori probabilities 
{tJ;}. Let denote T\ and A, two matrices such that Ei = 7}A,-7]- , A,- > and T{T^ = TLy E ., the 
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projection onto the support ofEj, for all 1 < i < n. Then necessary and sufficient conditions for a 
measurement {E^}, k = ?, 1, . . . ,n to be an optimal USD measurement are that there exists Z > 
such that 

ZE? = (2.15) 
EiiZ-rjipijEi = 0, l<i<n (2.16) 
T+iZ-riiPi)^ 1 > 0,1 <i<n (2.17) 

We could rephrase this theorem for two mixed states only. The statement then is slightly 
simpler. 

Theorem 4 Let po and pi be two density matrices with a priori probabilities r/o and 7]i. We 
denote by Pq~ and P^, the projectors onto the kernel of po and p\. Then necessary and sufficient 
conditions for an optimal measurement {E^}, k = 7,0,1 are that there exists Z > such that 



ZE? = 0, (2.18) 

E (Z-r]op )E = 0, (2.19) 

£i(z-i7ipi)£i = 0, (2.20) 

Pt(Z-T) Q p Q )Pt > 0, (2.21) 

/^(Z-mPi)^ 1 > (2.22) 



One could try to find the general solution for unambiguously discriminating two mixed states by 
solving the above conditions. However, in the general case it appears difficult to find a positive 
semi-definite operator Z fulfilling those conditions. Before ending this section, we can notice 
that 

Tr(Z) = P^ cess . (2.23) 

Indeed Eqn.(2.19) is equivalent to \/Eo(Z — f]ipi) \/Eo = 0. Its trace leads to Tr(Z£o) = 
T7oTr(p £ )- Similarly Eqn.(2.20) yields Tr(Z#i) = TjiTr(pi£i) so that Tr(Z£ ) + Tr(ZEi) = 
Psuccess- The completeness relation 1 = E? + Eq + E\ together with Eqn.(2.18) gives Tr(Z) = 
Psuccess- Later in this thesis, we will use Eldar's necessary and sufficient conditions to derive a 
theorem about the rank of the POVM elements of an optimal USD measurement and a new class 
of exact solutions of USD. 

2.5 Bounds on the failure probability 

2.5.1 Fidelity 

The fidelity F(po, pi) = Tr( yj i/popi y/po) is a quantity to distinguish two mixed quantum 
states po and pi. 



2.5. Bounds on the failure probability 
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We can consider the two extreme cases po = Pi and 5^ PQ _L y px . On one hand, if po = Pi then 
F(po, pi) = 1. On the other hand, if po and pi have orthogonal supports then F(pQ, pi) = 0. 
The fidelity takes value in [0, 1] . when F = 1, the two states are identical. When F — 0, the two 
states have orthogonal supports. It is not obvious that the fidelity is a symmetric quantity in its 
two arguments, though it is as we will show here ll40ll4Ti . We can first consider the fidelity of 
two pure states. 

F(|Y >(Y |, |Yi)(Ti|) = Tr( v /|Y )(Y || Y 1 )(Y 1 || Yo)(Y |) (2.24) 

= |(Y |Y 1 )|Tr( v /|Y 1 )(Yi|) 

= l(Yo|Yi>|. 

The fidelity of two pure states simply is the modulus of the overlap between those two pure 
states! The fidelity is here clearly symmetric. If we now consider mixed states, we can define 
the operators Fo = y^/PoPu/Po and F\ = J ^fp[pQ^fp{. They actually come from the polar 
decomposition 

y/fry/pi = F V = VF l . (2.25) 

As written in Eqn. (l2.25l) . the two operators Fq and F\ are unitary equivalent and their trace are 
equal. In other words, 

F{p l ,pj) = Tx{F l )=Tx{F ] ) (2.26) 

and the fidelity is symmetric. It might be sometimes difficult to work with the fidelity because 
of the three square roots involved in its definition and because of the noncommutativity of the 
density operators. For a review of its properties, the interested reader should look at Jozsa's 
1994 paper BUI inspired by Uhlmann's transition probability pT|. Let us however note here 
that in our work, the fidelity is given by F(pi, pj) = Tr ( ■J yfpipj y/pl) and not by F{pi, pj) = 
{Tr(-^/ y/plpj^/pi)} 2 11401 though the properties remain intact. 

Actually one can construct a distance measure from the fidelity, the Bures distance 
dg ures (pi, Pj) =2(1— F(pi, pj)). It is well know that the problem of minimum error discrimina- 
tion between two mixed states is linked to the trace distance as P err or = \ ( 1 — Tr( | T]opo — V lP 1 1 ) • 
As we are going to see through this thesis, a link between Fidelity and the failure probability Q in 
USD does exist. It is not as strong as the link between the trace distance as the error probability 
Permr in MED. In chapter 4, 5 and 6, we will intensively use the fidelity. 

2.5.2 Lower bound for the unambiguous discrimination of n mixed 
states 

Y. Feng et al. obtained a very general lower bound for unambiguously discriminating n mixed 
states {pi} with a priori probabilities {r],} [|42l . 
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Theorem 5 Let {p,} be n density matrices with their a priori probabilities 7],-. We define the 
fidelity of two states pi and pj as F [pi, pj) = Tr(-J \fpip \yfp~i)- Then, for any USD measurement 
a lower bound on the failure probability Q is 

2 > £mijF 2 (pi> Pj)- (2-27) 

V n &j 

Let us note here that another lower bound on the failure probability was derived by Y. Feng 
et al. (two of the three authors of Ref. 11421 ) in an unpublished work ll43l . Let us notice that this 
bound is given as an upper bound on the success probability. 

Theorem 6 Let {p, } be n density matrices with their a priori probabilities {fy}. First we define 
the subspace Mix(pi) as Mix(pi) = S^p i ^^Lj^i^pj- Second, we divide each pi in two parts, pi 
and pi such that 5^'p i = Mix(pi) and 5?p. fl S^p i = 0. Finally we define the fidelity of two states pi 
and pj as F(pi, pj) = Tr{^J ' y/plPjy/pi)- Then, for any USD measurement an upper bound on 
the success probability P success is 

Psuccess < £ iMpi) ~ . / —3- £ mVjF 2 (fa Pj) ■ (2-28) 

i=i V ¥J 

This last bound is tighter than the one in Theorem 5 since YJi=\ r lf^ T {pi) ^ 1- The equality holds 
only if the density matrices p ; do not have common subspaces. In that case, the two lower bounds 
in Eqn.(2.27) and Eqn.(2.28) are equal. We now focus on USD of two density matrices only. 
Rudolph et al. derived both lower and upper bounds on the failure probability to unambiguously 
discriminate two mixed states. This is the object of the last subsection of this chapter. 

2.5.3 Lower and upper bounds on the failure probability for the un- 
ambiguous discrimination of two mixed states 

Lower bound 

In Ref. 11261 . Rudolph et al. derived their lower bounds considering some purification of the two 
mixed states po and p\. Moreover, an interesting property of the fidelity is the following. Given 
two mixed states, we can consider all their possible purification and their overlap. In fact, the 
fidelity equals the maximum of the modulus of those overlaps. It is therefore not surprising that 
those lower bounds involve the optimal failure probability of two pure states where the overlap 
is replaced by the Fidelity (see Fig. 12.31) . More precisely, we end up with 

Theorem 7 Let po and p\ be two density matrices with a priori probabilities r/o and r\\. Let 
define the fidelity F = Tr( y^/pop~u/po) between these two mixed states. Then a lower bound 
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on the failure probability of unambiguously discriminating po and p\ is 



2° pt > m + TfoF 2 for \ —<F, 

V 

2 opt > 2yrPjTFforF< A /^<4, 

V ^0 r 

Q opt > rio + riiF 2 for ^< 




(2.29) 
(2.30) 
(2.31) 







Qopi > m + mF 





Q° vt > no + '). F 



Figure 2.3: Lower bounds on the optimal failure probability for USD of two density matrices 



Upper bound 

In the same paper ll2~6ll . the authors presented an upper bound on the optimal failure probability 
for unambiguous discrimination of two mixed states. This bound comes from considering several 
two dimensional USD problems rather that a global USD problem. The eigenbases for Eq and E\ 
here depend only on the supports of po and p\ and not on their eigenvalues. This leads naturally 
to an upper bound on the failure probability since the eigenvalues of po and p\ would allow to 
refine the measurement. The theorem presents a lower bound on the success probability instead 
of an upper bound on the failure probability. 

Theorem 8 Let po and p\ be two density matrices with a priori probabilities r/o and r\\. We 
denote the dimension of their kernel and J??i by sq and s\ and assume that sq > s\. There 
exist orthonormal bases {\k J h )}J =l for (b=0,l) such that for 1 < j <sq, 1 < i < s\, 



(k J \k\)=Cos(e j )8 ij , 



(2.32) 



where the 0j are the canonical angles between J#6 and ,J^\. In this case, a lower bound on the 
optimal success probability Process is 



so 



^success — $2 ^success ( I ^0 ) ' I M ) ) 12 (^0 1 Pi I ^o) ' 



(2.33) 
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2. Optimal Unambiguous State Discrimination 



where 

«» = { ^i- 2C -W V^i &r C«(9;) < J± (234) 
^ A^ox otherwise 

wif/i Aq = Tj (*i|po|*j), A{ = T]i(^|pi|^), A/ n . 3 = ra/n{A^Aj} andA J max = max{A J Q ,A[}. 

Let us note that we will detail the construction of such orthogonal bases {\k J b )}J = i in Chapter 3 
when we will present the optimal unambiguous discrimination of two subspaces. 

In the next chapter, we will find that any USD problem can be reduced to some standard 
situation. We will then see that some important tasks in Quantum Information Theory which are 
related to the USD of some mixed states can actually be reduced to some pure state case. 



Chapter 3 

A standard form 



We are searching for an optimal USD measurement to discriminate two arbitrary density 
matrices po and p\ with a priori probability T}q and rji respectively. We find that this general 
problem can be reduced to a simpler standard situation thanks to three reduction theorems 
dealing with simple geometrical considerations. As their names indicate, the three reduction 
theorems allow to reduce the dimension of the USD problem. In fact, the reduction can also be 
applied to the case of more than two density matrices. 

It is important to notice here that all the results on USD of mixed states known so far 
are reducible to some pure state scenarios. These cases are state filtering, unambiguous 
discrimination of two subspaces and unambiguous comparison of two pure states. Those three 
cases of USD of mixed states can be solved using some reduction theorem and the result of 
Jaeger and Shimony about USD of two pure states only. This underlines the fact that those cases 
were solved first because no new techniques were needed. In the following we will often refer 
to non-reducible mixed state case as generic USD problem. In the next chapters we are going 
to present two classes of exact solutions for such generic USD problems. But first of all, let us 
present, prove and use the three reduction theorems. 

The first reduction theorem states that, if two density matrices share a common subspace 
(see Fig. 13.11) . no unambiguous discrimination is possible on it. Indeed any state vector in such a 
common subspace belongs to both po and p\ so that no conclusive result is possible. The failure 
probability restricted to this common subspace then equals unity. There is no optimization to 
perform onto this common subspace and we can focus our attention on the USD problem onto 
the orthogonal complement of this common subspace. 

The second theorem is easy to understand, though the proof happens to be subtle. Let us 
consider the support S^ Po and S fi Pl of two density matrices. Let us assume that there exists a 
subspace of 5? Px orthogonal to y po (see Fig. 13 .21) . This subspace can be equivalently denoted by 
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3. A standard form 



J7 pi Pi S^p^ or J/p, fl ^p - If we perform any measurement on that subspace, we can only detect 
pi but never po since the measurement is orthogonal to y po . The difficulty step is to see that 
such a strategy is optimal. Here again, no optimization onto the subspace y pi fl J£p is needed. 
After splitting off y pi fl Jf Pa , we are left with a smaller USD problem. Of course, a similar 
reduction can be performed for the subspace y po fl Jt pi . 

The last theorem refers to some block diagonal structure of the supports y po and 5? Px of our 
two density matrices po and p\. If the supports y p(i and y pi can be simultaneously decomposed 
into a direct sum of some subspaces, it seems reasonable that the optimal measurement can have 
the same property. Moreover we can choose the optimal measurement onto the total Hilbert 
space to be the direct sum of optimal measurements onto the smaller subspaces. In other words, 
we only have to look for optimality on each orthogonal subspace. This again simplifies the 
optimization task. 

Let us now derive the three theorems. 



In the first theorem, we will consider the situation where the supports of the two density matrices 
have a common subspace. This is the case whenever we find that 



Here M 3 is the Hilbert space spanned by the two supports. In this case, it can be written as 



where J^n = yp ^ 1S me common subspace of the two supports, and J^", its orthogonal 
complement in Jt? (see Fig. 13.11) . The first reduction theorem will eliminate the common 
subspace J^n from the problem. The intuitive reason is that in this subspace no unambiguous 
discrimination is possible, so the population of the two density matrices on it will contribute 
always only to the failure probability, never to the conclusive results. This is made precise in the 
following theorem. 



3.1 Overlapping supports 




(3.1) 



(3.2) 



3.1. Overlapping supports 
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Figure 3.1: Illustration of a common subspace between po and p\ 



Theorem 9 Reduction Theorem for a Common Subspace 

Suppose we are given two density matrices po and pi in Jtf? with a priori probabilities t\q 
and r\\ such that their respective supports S^ Po and 5^ Px have a non-empty common subspace 
J£n. We denote by M" the orthogonal complement of in while Hje n an d n 
denote respectively the projector onto and M" . Then the optimal USD measurement is 
characterized by POVM elements of the form 

E opt = E 'opt (33) 
E opt = E 'opt (34) 

ET = E^ + n^ (3.5) 



where the operators E^ pl ,E^ pl ,E 7 0pt form a POVM {E^ pt } with support on <!%" describing 
the OptUSDM of a reduced problem defined by 

p6 = ^n^p n jr , t^=» No = Tt(poUj ei ) (3.6) 

p; = i-n^ Pl n^, v[ = ^, m = Tr^n^) (3.7) 

N = N rio + N l ri l . (3.8) 

And finally, the optimal failure probability Q opt can be written in terms ofQ' opt , the optimal 
failure probability of the reduced problem, as 



Q°V = i_Af + jV2 /opt . 



(3.9) 
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3. A standard form 



Proof To prove the reduction theorem, we first need to recall that a USD measurement de- 
scribed by the POVM {E{\ satisfies Tr(£opi) = an d Tr(.E'ipo) = by definition. It means, as 
a consequence of Lemma 1 given in the previous chapter, that Se -L Sp, and 5^, J_ S Po . Since 
J#h is a subspace of 5 Po and S Pl , it follows that Se -L and Se { -L J#h- Therefore, by writing 
the block-matrices in Jt? = J% © M" , we have 





E' Q 



E[ 



(3.10) 
(3.11) 



The completeness relation on M' implies firstly 

X t E' 7 ) =Uj « +E l (3 - 12) 
and secondly by the completeness relation on the reduced subspace M" 

= (3-13) 

k 

It follows also that the operators E' k (k = 0, 1,?) are positive semi-definite operators. Therefore, 
by definition, {E'A is a POVM on M?' . The fact that E-> is equal to identity in the subspace J^n 
is here a direct consequence of the property of an USDM on ffl . Next we will see that {E' k } is a 
POVM of a USD in M" . 



We define Tl r ^ n and Yl^t as the projector onto J^n an d M" respectively. Thus FIj^ © 
Yljfi = \jf. For any USDM, because of the diagonal block form of the POVM, we find for Q 

Q = 77oTr(p £ ? ) + 77iTr(pi£ ? ) 

= (l-iVo)T?o+(l-M)T?i (3.14) 
+ (iV T] +iV 1 T7i)(T7oTr(p^) + ^Tr(^^)) 

wimp/, = -—, — — — r-TI poTLjif" (3.15) 
P0 Tr(p n^) ^ H ^ 

Pi = t ( n It \ n ^'Pl U ^'- ( 3 - 16 ) 

Here r/- (z = 0, 1) is the a priori probability corresponding to the new density matrix p- (t]q + 

i! = D 

^ = 7t£^- N ° = Tl(p ° n ^ <3 ' 17) 
* = d™'" 1 = Tr(pirw ' ai8) 



3.2. Trivial orthogonal subspaces of the supports 
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We notice that y p i D y p i = 0. Moreover, Tr(p £i) = implies Tx{p' Q E[) = and 
Tr(pi£ ) = implies Tr(p[E' Q ) = 0. Then {E' k } defines a POVM describing a USDM on 
{p[, rj-} in M" . The problem is now reduced to the subspace M" . We now focus our attention 
on the optimality of the reduced USDM. 

We can write Q as 

Q = (l-iVo)T7o + (l-iV 1 )r]l + (^oT]o + iV 1 T] 1 )e , (3.19) 
= l-N + NQ' 

where Q' = t]qTi:(PqE' 7 ) + r}[Tr(p[E' 7 ) is, by definition, the failure probability of discriminating 
unambiguously p' Q and p[ in 3%" with a priori probabilities 7]q, r)[. 

The previous equality implies that the failure probability Q is minimal if and only if the 
failure probability Q' is minimal. Thus we have that {E^} describes an optimal USDM on 
{pi, T]i} <^ Q is minimal <^ Q' is minimal <=> {E' k } describes an optimal USDM on {p/, rj-}. 
This completes the proof. ■ 

Let us note here that two subspaces that do not have a common subspace are not necessarily 
orthogonal. The formal statement is y po D y pi = {0} y po _L y pi . Moreover we can give an 
easy way to know whether the two supports overlap of po and pi. In fact, it suffices to check 
whether the equation dim(J$?) = rank(po) + rank(pi) = rank(po + pi) holds. Marsaglia and 
Styan proved that additivity of rank of two matrices is related to the intersection of their column 
and row spaces in a simple way ll44ll . Their result is given in the following theorem. 

Theorem 10 Let A and B be two complex mxn matrices. Let % and ^ be their column spaces 
and Ma and Mb, their row spaces then 

rank(A + B) = rank(A) + rank(B) if and only if dim^A^^B) = dim{&A^&B) = {0}. 

In the more restricted case of two density matrices, which are Hermitian matrices, the column 
and row spaces simply are the support ^ p = M p = y p . 

3.2 Trivial orthogonal subspaces of the supports 

We now consider the case where the supports of the two density matrices have no common 
subspace. That can always be achieved thanks to the previous reduction theorem for common 
subspace. If there is a part of y pi orthogonal to yp , we can decompose y pi into this subspace 
and another one (see Fig. 13.21) . It turns out that this subspace of y pi orthogonal to y po can be 
split off and leads to an unambiguous discrimination without error. The same is true for y po . 
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3. A standard form 



Theorem 11 Reduction Theorem for Orthogonal Subspaces 

Suppose we are given two density matrices po and p\ in with a priori probabilities t\q 
and r\\. Assuming that their supports y pQ and y px have no common subspace, one can 
construct a decomposition 

M' = M 3 '®M 3 ' L (3.20) 

with jtf^ = s^@ s\, = ^p n and = n y Po . 

The solution of the optimal USDM problem can be given, with help of JJ „>± and II the 
projection onto and Sf^, respectively, in = M" © ffl , by 

E / = E'° pt + n K (3.21) 

E opt = E 'opt +u (322) 
E opt = E 'opt_ (3 23) 

The operators Eq P ' ,E° pt ,E° pt form a POVM {E k ° pt } with support on M" describing the 
OptUSDM of a reduced problem defined by 

p' Q = ^-U^p Q n^, 77^= « V = Tr(p n^) (3.24) 

p[ = ±-n^, Pl n^, t?; = ^, M = Tr( Pl n^) (3.25) 

N = N r] + N l r ll . (3.26) 

And finally, the optimal failure probability Q opt can be written in terms ofQ' opt , the optimal 
failure probability of the reduced problem as 

gopt = N Q'0 P t_ (3 _ 27) 



Proof We translate the problem using a Naimark extension and a projection-valued measure 
(PVM). This idea is inspired by the first work of Sun et al. Il32l where an extended Hilbert space 
has been used. Let us repeat the Naimark theorem. 

Given a POVM {E k } on a Hilbert space J4f, it exists an embedding of Jt? into a larger Hilbert 
space M such that the measurement can be described by projections onto orthogonal subspaces 
in More precisely, there exist a Hilbert space an embedding $ such that $J4? = and a 
PVM {R k } in M such that with P, the projection defined by P& = Jf,E k = PR k P, Vfc. 



3.2. Trivial orthogonal subspaces of the supports 
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To the three POVM elements in M* correspond three PVM elements Rj, in The Hilbert 
space ^ can be decomposed into orthogonal subspaces 

^ = S%®S fi Rl ®S f R 1 (3.28) 

which give raise to non-orthogonal subspaces in Jf? as S^E k = P^R k P- We can therefore translate 
properties of the USD POVM to the embedding of J? into St. 

Next we take a look at the embedding of ,5^ P(S and S fi Pl into M and we translate the conditions 
for an USDM into the embedded language. We denote the embedded subspaces of M by the 
same symbol as the original subspace of M '. We can here introduce the projector P 1 - onto the 
orthogonal complement Jf? 1 - ofJ^mM(P + P ± = Since y po G Jtf, we have Tr(p i?i ) = 
Tr(PpoPRi) = Tr(po£i) = 0. This implies that y pi) is orthogonal to S fi R v Similarly, we find 
that y pi is orthogonal to =5^ . Therefore, we can write 

^p CJ^ ©J^ ?0 (3.29) 

y Pl cy Rl ®y Rn (3.30) 

where J^ ?0 and S fi R n are defined as subspaces of 5^r 7 with minimal dimension fulfilling the 
above decompositions in the sense that S^r 1a = Support(n y R y Pi TI ,y Ri ) for i — 0, 1. 

The optimality condition means in particular that no information should be obtained from the 
conditional states following an inconclusive result. If the two failure spaces y R?0 and 5?R n are 
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3. A standard form 



different, it will be possible to distinguish the conditional states which arise from a projection 
onto y^ ll32l . Indeed a detection in an orthogonal direction to one of the two subspaces will 
tell us which failure space was it or equivalently which state was sent. Therefore the optimality 
condition implies that ^? ?0 = <^R n and then 

= ^r ?0 = ^R? r (3.31) 

This is an important necessary condition for the optimality of a USD POVM. In the framework 
of the Naimark extension, this condition translates as follows. The equality of J?/f ?0 and ^R n 
implies that a subspace y^~ = Jf Po fl y pi satisfies y§ C S^r x in order to assure that the overlap 
between any state in y^ and any state in y po will be zero. Similarly, y^ = J% Vl fl y po C <5^ . 
Then there exist two subspaces in y& x and ^f> in y^ such that 

y Rl = y^®-m (3.32) 
y Ro = yt®<%- (3.33) 



The orthogonal projection R\ then can be decomposed into a sum of orthogonal projectors 
as ITyx + lljfj, with Ily±H^ = 0, and the orthogonal projection Rq as H^± + TLj%, with 
Ily±Tlj^ Q = 0. These projectors are mapped into J$? via the projection P. Since y^~ is already 



in ffl, we have PW^^P 



n^x. We define E\ 



PUj^P, Vi = 0, 1 so that 



E' + T1 

E[ + n 



y, 1 



(3.34) 
(3.35) 



Furthermore, the two supports y E i and y^ are orthogonal since Tl^ Q Tly± = implies 
U^ PPYly±P = so that PYl^PPYly^P = E' TIy± = 0. Similarly the two supports y^ 
and y^~ are orthogonal too. 

Moreover, y E() ± y pi and y§ G y pi so that y Eo _L y^ . Similarly, we have y Ei _L y^. 
Then E' Q and £| have support on a subspace which is the complementary orthogonal sub- 



space of je 3± = yj-®y 1 A 

Therefore in = JT' ( 



1 ®y^ = ^' l ®^' A -, we find 







V 



o 







1 







1 





\ 



(3.36) 



(3.37) 



3.2. Trivial orthogonal subspaces of the supports 
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From here, we will follow the same argumentation as we used in the proof of Theorem 9. The 
completeness relation on ffl implies firstly 

£^0 



and secondly the completeness relation on the reduced subspace M" 

E4 = 1jt. (3-39) 

k 

It follows also that the Et (k = 0, 1,?) are positive semi-definite operators. Therefore, by defini- 
tion, {E' k } is a POVM on M" . 

Let us note that y po c © © y R?Q and y pi c ytf~ © M{ © y Rn . The fact that yj- c 
,5^ implies that 

^Po = ^©^ (3-40) 
with Jf%cJ#o® y R?0 . In the same way, with J^/ C^i® y Rn , 

«^ Pl =^j L ©X- (3-41) 
Therefore, we can introduce a reduced problem onto defined such that 3? = Jt" © y§ © 

y x \ 

For any USDM, because of the diagonal block form of the POVM, we find for Q 

Q = 77oTr(p £ ? ) + 77iTr(pi£ ? ) (3.42) 

= (iV T] +iV 1 T7i)KTr(p^) + T] 1 Tr(p;^)) 

with Po = w— i=f — rll^poU^ (3.43) 
Tr(p n^) 

^ = T/,n ''//'Pi"//'- (3-44) 
Here T]- (z = 0,1) is the a priori probability corresponding to the new density matrix p[ 

(io+i! = D 

^ = d™' A ' = Tr(A>IW <345) 

Moreover, Tr(poEi) = implies Tt(PqJ?j) = and Tr(piEo) = implies Tr^i?^) = 0. Then 
{E' k } defines a POVM describing a USDM on {p/} in M" . 
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We can rewrite the failure probability Q as 



Q=(N Tlo + NiTii)& (3.47) 



where Q' = 7]qTt(PqE' ? ) + J][Tr(p[E' 7 ) is, by definition, the failure probability of discriminating 
unambiguously p' Q and p[ in Jrff' with a priori probabilities r\' Q and r\[, respectively. 

And again, we have that {E^} describes an optimal USDM on {pi, 7],} <^> Q is minimal <^ 
Q' is minimal <^> {E' k } describes an optimal USDM on {p-, r\[}. This completes the proof. ■ 



3.3 Block diagonal structure 



It is possible to state a last geometrical theorem which deals with two block diagonal density 
matrices po and pi. Schematically, po and pi are then of the form 



/□ \ 

o □ o 
V o o □/ 



The problem of unambiguously discriminating such two density matrices can be reduced to 
smaller USD problems onto each one of the orthogonal subspaces. This is made more precise in 
the next theorem. 



3.3. Block diagonal structure 
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Theorem 12 Reduction Theorem for two block diagonal density matrices 
Suppose we are given two density matrices po and p\ in with a priori probabilities 
t\q and r\\. Suppose that po and pi are block diagonal (in other words, it exists a set of 
orthogonal projectors {Tl k } such that Y! k= \ = 1 ond pi = Y! k= \ ^kPi^h i = 0, 1 ). Then 
the optimal USD measurement can be chosen block diagonal where each block is optimal 
onto its restricted subspace. 

More precisely, the optimal USD measurement is characterized by POVM elements of the 
form 

E T = Y, E t opt - < 3 - 48 ) 

k 

For k = l,...,n, the operators Eq opt ,E kopt ,E k opt form a POVM {E kopt } with support on 
Sfp k describing the OptUSDM of the reduced problem defined by 

p k =^ll k p Q U k , i?J = ^r. 4 = Tr(p n,) (3.49) 

P f = ^n, Pl n,, ??f = ^, N k = ix{ Pl u k ) (3.50) 

N k = Nfa +N$T]i . (3.51) 

And finally, the optimal failure probability can be written in terms of Q k opt , the failure 
probability of the reduced problems, as 

e opt = "L N kQT- 0.52) 

k 



Proof We start with two block diagonal mixed states po and pi with a priori probabilities 
r/o and r\\. In other words, we assume that it exists a set of orthogonal projectors {IT*} such 
that Y!l=\ = 1 and p,- = ££ =1 Yi kP jYl k , i — 0, 1. Next, we denote <9xi k , the support of the 
projector FL;. We first show that only the restriction of the POVM to the n orthogonal subspaces 
S^u k is relevant to the failure probability. Then we will show that optimality on each orthogonal 
subspace S^iu leads to optimality on the total Hilbert space. Let us consider a USD POVM {Ej} 
onto Jf? and its failure probability Q which can be written 
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Q = J>Tr(£ ?P( ) (3.53) 

i 

= £T7/rr(£ ? (£n, P/ n,)) 

i k 

= £J>Tr(n fe £ ? n fe n fePi -nfc) 



k i 

We can obviously define n reduced density matrices onto the n subspaces ^S^fr. as 



rf - ^ 0.54) 



/Vf = TrCn^Pi) (3.56) 

with N k = Y.i^ili- We can also consider the restrictions of the POVM elements Eq,E\ and 
E? onto those n subspaces. Thus 



E k = n k E YI k (3.57) 



El = U k E x Il k 

e\ = n k E ? n k . 



Obviously those operators Ef (i = 0, 1,?) are positive semi-definite and add up to YI k since 
Y.iEi = 1. Each restriction onto ^\\ k of a POVM \E\\ then forms a POVM onto the subspace 
y Uk . Moreover Tr(£*p*) = Tr(U k E lP] Il k ) = Tr (11^11*)$, since E iPj = E lPl 8,j for i,j = 
0, 1, so that the n POVMs are n USD POVMs. 

As a consequence, the failure probability for any two block diagonal density matrices can be 
expressed in terms of the failure probabilities Q k = £ ; r\ k Tr(E k pf ) of the n reduced problems as 

Q = Y, N kQ k - (3-58) 

k 

We can now show that if each block is optimal then the block diagonal POVM onto Jt? is optimal 
too. 

To prove it, let us consider an optimal USD POVM onto each one of the n orthogonal sub- 
spaces <¥n k - We denote Q ko P t the optimal failure probability onto ^n k - By definition of the 
optimal failure probability, Q k > Q ko P' for each subspace J?n k - Since both N k and Q k are posi- 
tive numbers, this yields 

Q>Z N kQ k ° Pt - ( 3 - 59 ) 

k 



3.4. A standard form of USD problem 
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This bounds can be reached for {Ej} being the direct sum of the n optimal USD POVMs {Eft 
i.e. Ej = Y<k=lEj> J = 0, 1, ?. The completes the proof. ■ 



3.4 



A standard form of USD problem 



At this point, it is useful to introduce a notation to summarize our knowledge about the USD 
of two density matrices. We have = y po + y pi then dim (Jt?) = dim (y po ) + dim (y pi ) — 
dim(^p ny pi ). It implies, by denoting the dimension of the Hilbert space Jt? as d, that the 
respective ranks tq and r\ of the density matrices po and p\ satisfy 



For example, the case of two density matrices of the same rank (n — 1) in an Hilbert space of 
dimension n described by Rudolph et at ll26l can be written as "(n — 1) + (« — 1) > n" while 
the USD between one pure state and a mixed state described by Bergou et al. 11321 133ll3"4ll can be 
characterized as the "1 + n = (n + 1)" case. We will see in the following section that important 
tasks in quantum information theory can be solved elegantly thanks to those three reduction 
theorems. 

First of all, let us discuss some immediate consequences of the three above theorems. The 
first reduction theorem corresponds to the elimination of the common subspace. A common 
subspace is present when tq + r\ > d holds. Its dimension is d n = r$ + r\ — d. Therefore, after 
elimination of that subspace, we end up in the case r' + r[ = d' with r' Q = tq — J n and similarly 
for r[ and d' . Then, we can reduce the Rudolph's case of discriminating unambiguously two 
density matrices of the same rank (n — 1) in an Hilbert space of dimension n to the "1 + 1 = 2" 
case of two pure states because the common subspace is (n — 2)-dimensional. Rudolph et al. Il26l 
already noticed it in their paper. The reduction is constructive given po and p\. 

The second reduction theorem corresponds to the elimination of the orthogonal part of one 
support with respect to the other, i.e., J^ Pl^p, and Jif Px r\y po . The non-empty subspaces 
J^ Po fl y pi and J^>, fl y po can be found systematically. For example, Jtf Po PI y pi can be found 
by projecting y po onto y pi and then by taking the complementary orthogonal subspace in y pi 
of that projection. As a matter of fact, this assures that we can reduce a general USD problem 
always to that of two density matrices of the same rank r, r < min(ro,ri), in a Hilbert space of 2r 
dimensions. Indeed, if after the first reduction, the rank of p[ is bigger than the rank of p' Q , then 
the subspace J^ p i fl y p i is at least of dimension r[ — r' and can be eliminated. With the help of the 
first two reduction theorems, we can reduce any problem of discriminating unambiguously two 
density matrices po and pi, with rank tq and r\ respectively, in a Hilbert space Jrf? , into a problem 
of discriminating unambiguously two density matrices p' and p[ with rank r (r < min(ro,ri)) 



ro + n > d. 



(3.60) 
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3. A standard form 



in M" C 3^, a 2r-dimensional Hilbert space. The reduction is constructive. The first theorem 
allows us to split off the common subspace and the second theorem leads to the reduce problem 
of discriminating unambiguously two density matrices of the same rank. The third theorem tells 
us that if the two density matrices have a block diagonal structure, we can reduce the problem of 
unambiguously discriminating them to some smaller ones, each one corresponding to a block. In 
fact, the three reduction theorems allow us to define a standard form of USD problem as follows. 

Definition 5 Standard form 

Any Unambiguous State Discrimination problem of two density matrices of rank ro and 
r\ is reducible to that of two density matrices of the same rank r < min(ro,n)J in & 2r- 
dimensional Hilbert space without overlapping supports, without trivial orthogonal sub- 
spaces and without block diagonal form. Such a problem is called a standard Unambigu- 
ous State Discrimination problem. 

The expression 'trivial orthogonal subspaces' stands for the subspaces J€ P{) ny pi and Jt Pl fl 
y po . It is also interesting to note that the dimension of the failure space can not be greater 
than the lowest rank of the involved density matrices. In the language used in the proof of the 
second reduction theorem, we first have E? = PR?P so that dim( t 5^£ ? ) < dim(=5^ ? ). Second the 
dimension of yR 7i can not be greater than r, because =5^,. = support^?^,^?), for i = 0, 1, and 
= = yR ?r Therefore dim=5^ ? < min/dim^p, and we can define the maximum rank 
of E? as 

r'^ = min{r Q , ri ). (3.61) 

This result looks natural considering that we can finally reduce any problem of discriminat- 
ing two density matrices with rank ro and r\, respectively, to the problem of discriminating two 
density matrices of the same rank r, r < min,r/. 



Finally, a generalization to more than two density matrices can be achieved. Considering 
n density matrices p^ (k — 0...n — 1) with a priori probabilities 7]^, we can construct n pairs of 
density matrices 

po = pi, ie [0,..,n-l] (3.62) 

and 

= r; V7/P/ 

with f)o = T]i, f)i = 1 — rji, and apply the two reduction theorems to these two density matrices 
in the following sense (notice that pi has no physical meaning). As soon as a common subspace 
between any y po and y pi exists, we can split it off from all the ypfs, because if we cannot 



3.5. Applications of the reduction theorems 
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discriminate unambiguously this part of the support of po and p\ then we can not discriminate 
unambiguously between this part of the support of all the p 7 . The second theorem must be used 
more carefully. As soon as a subspace of =5^p is orthogonal to 5?p x (J£p, PI ^p / {0}), we 
can eliminate it from the problem because it is orthogonal to the supports of all the pj, j ^ i. 
However we cannot eliminate a subspace of S?p x orthogonal to S e 'p (J^s PI 5?p x ^ {0}) because 
we know nothing about the orthogonality of this subspace for all the states in p\ . In other words, 
we can only reduce the density matrix p,- corresponding to po. 

In the following section we are going to apply the reduction theorems to three important 
tasks in quantum information theory. Those tasks are State Filtering, Unambiguous Comparison 
of two subspaces and Unambiguous State Comparison of two pure states. We are going to see 
that those three tasks are reducible to some pure state case only. 



3.5 Applications of the reduction theorems 

3.5.1 State Filtering 

Let us consider n pure states with a priori probabilities {pi}, i = 0, ...,n — 1. We may want 

to group them in several sets and to unambiguously discriminate among these sets. This task 
is called State Filtering Il32ll34ll . The simplest case deals with two sets only where the first set 
contains one pure state and the second set regroups the remaining n — 1 states. This problem was 
studied in various papers by Bergou et al. J32l [331 EH who gave the complete solution in ll34l . 
We derive here this last result is an extremely simple way thanks to the second reduction theorem. 

We have to unambiguously discriminate the two sets {|Yo)} and {|Y,)} ; =i,...,; 7 -i. We can 
consider the density matrices corresponding to these two sets as well as their a priori probabili- 
ties. The first density matrix obviously is po = |Yo)(Yo| with a priori probability t\q = pq. The 
second mixed state can be written as 

n-1 

Pi=J>/|Yi><Y/|. (3.64) 

This is not a proper density matrix since it is not normalized. We then must write p\ = 
' „'-i • Its a priori probability simply is r\\ = Y^=\ Pi = 1 — Po- State filtering finally 
is equivalent to unambiguously discriminate 

Po = |Y )(Y | (3.65) 

with a priori probability r/o = po and 
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with a priori probability r\\ = YJjZi Pi- 
After writing these two density matrices, the solution to the problem is trivial. 

Indeed a consequence of Theorem 1 1 is that we can reduce the problem of USD between 
a pure state and a density matrix, a "1 + n = (n + 1)" case, to the problem of discriminating 
unambiguously two pure states, a "1 + 1 = 2" case, by splitting off Jf Po fl y pi of dimension 
(n — 1 ) . The two reduced states are the original pure state | Yo) and the unit vector corresponding 
to the projection of po onto the support of the mixed state pi. This unnormalized vector is given 
by = rii|Yo), where III is the projector onto the support of p\. The corresponding unit 

vector simply is |Yg) ~ 



Ml 



Theorem 1 1 tells us that the optimal failure probability Q° pt for State Filtering is given by 

G°*=W*(|Yo),|Y<))), (3-67) 



with 



p' Q = ^-U^p Q U^, 77o = ^r> ^o = Tr(p n jr ) (3.68) 

p| = i-n^ Pl n^, TJi = ^. M = Tr( Pl n^) (3.69) 

N = N r l0 + N l r] u (3.70) 

jr' = {|Y ),|Yg>}. (3.71) 



Furthermore, the optimal failure probability for two pure states |Yo) and |Yq) with a priori 
probabilities r\' Q and r\[ is given by 



£ opf (|Y ),|Yo)) = + ^|(Y |Y^)| 2 for < |(Y |Yg)|, 

V "o 

C^CW.TO) =2 v /^|(Yo|Y' / )| if |(Y |Yo)| < J| < |(Yo | Yg)| , 



e° , ' f (|Yo),TO) = ^+T ] ;KY TOi 2 if 



A. 



KYoTOi - V *i6 



(3.72) 
(3.73) 
(3.74) 



therefore the optimal failure probability Q° pt of the non-reduced problem becomes 
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e^A^+T^KYoTOl 2 ) for W^<|(YoTO|, 

V 

=iV(2 v /^|(YoTO|) if |<YoTO| < M < J^f^y 



<2 opt = ^K+^i(Y TOi 2 ) if 



(3.75) 
(3.76) 
(3.77) 



If we denote S = £J~{ 1 (Y 1 Y,-) | 2 , we 



find 



iVo 



Hi 

|{YoTO| 



ml mi 


2 


rioNo _ 




N 


N 


7]i/Yi _ 


S 


N 


N\\¥^\\ 2 


1 1 W'l 1 

ll T oll- 





(3.78) 
(3.79) 

(3.80) 
(3.81) 

(3.82) 



We finally end up with 



gopt 

gopt 
gopt 



po||Y^|| 2 + ^^if _^<p , 

I lw//| |2 1 1 ait// 1 14 

ll T oll ll T oll 



2^f~PQ\fS if 5 < po < 
Po + S if po < S- 



|w//| |4 ' 

l T oll 



(3.83) 

(3.84) 
(3.85) 



3.5.2 Unambiguous Subspace Discrimination 

To unambiguously discriminate two subspaces, one has to unambiguously discriminate their 
respective bases. We can therefore consider the two ensembles corresponding to these two bases 
with a flat distribution because the basis vectors all possess the same probability of appearance. 
In fact we consider the projectors onto those respective bases as unnormalized mixed states and 
try to unambiguously discriminate them. In that sense, subspace discrimination is a special case 
of mixed state discrimination where the two density matrices are proportional to the projectors 
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onto the respective subspaces. 

There is a infinite amount of basis in which one can write a projector. Therefore the difficulty 
is to find a suitable basis of the space spanned by the two subspaces to discriminate. Such 
a suitable basis is given by the so-called canonical bases which allow us to write the two 
projectors in a block diagonal form, where each block is two-dimensional. This technique was 
used by Rudolph et al. for the derivation of the upper bound on the failure probability Q. Thus 
the unambiguous discrimination of two subspaces can be reduced to some pure state case and, 
because of that, be solved. 

First, let us repeat that the first two reduction theorems permit us to focus our attention on the 
unambiguous discrimination of two subspaces So and Si of rank r in a 2r-dimensional Hilbert 
space. Next we choose an orthogonal basis {|a,-} } of So and an orthogonal basis {\bj) } of Si . The 
unambiguous discrimination between these two subspaces then corresponds to the unambiguous 
discrimination of po = }L/ and pi = \bj)(bj\. 

Given two subspaces So and Si, it is always possible to find an orthonormal basis of 
So and an orthonormal basis {\bj}} of Si, called canonical or principal bases such that (ai\bj) = 
Cos(6j)8ij, Cos(Oj) > 0. In such a basis, the projectors onto So and Si are decomposed into a 
direct sum of r two-dimensional subspaces. Thanks to theorem 12, the optimal solution to USD 
of two pure states is the only requirement for an optimal unambiguous discrimination of So and 



In fact, we can assume without loss of generality that (ai\bj) = Cos(0j)8jj, Cos(Oi) > 0. In- 
deed, we can always construct the so-called canonical bases and {\bj)} for two subspaces 
if we follow Rudolph's technique ll26ll . LetXj. be the (2r)xr matrix whose columns span S^. We 
then write a singular value decomposition of xlX\, 



where the t4's are two rxr unitaries and S is positive semi-definite and diagonal with Su = 
Cos(6i), (6 E [0,2n]). Let us define the vectors \ai) as the i column of XqUo and the vectors 
\bj), the / column of X\U\. The set respectively {\bj)}, forms an orthonormal basis of 

So, respectively Si, since it is merely a rotation of a former basis. Moreover the vectors |a;) and 
\bj) satisfy (aj\bj) = Cos(0j)8jj. The angles 0j are called the canonical angles and, the vectors 
\ai) and \bi), the canonical vectors. \aj) and \bt) together span the total Hilbert space. The funda- 
mental property (ai\bj) = Cos(Oj) 8jj allows us to write po and pi in a block diagonal form, where 
each block is spanned by {|a ; ), \bj)}. Indeed, in the basis {\a\) \b\), \a,2), \bi),- ■ -,\a r ), \b r )}, the 
two density matrices po and pi takes the form 



Si. 




Xi = u su 



•t 



(3.86) 
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/□ \ 

p k = o n ° 

V o o □/ 

where, each block is a two-dimension subspace spanned by {|a;) \bi}}, orthogonal to the r— 1 
other two-dimensional subspaces {\a k ) \b k )}, k—\,...,i—\,i-\-\,...,n. 



Thanks to theorem 12 we can express the failure probability of unambiguously discriminating 
Sq and Si as 



gopt _ ^jykQkopt^ 



(3.87) 



where the Q ko P l are the optimal failure probabilities for unambiguously discriminating \a k ) and 
\b k ) with their corresponding a priori probabilities and r\\. 

We can easily calculate all those quantities where Tl k is the projector onto the two dimen- 
sional subspace spanned by \a^) and \b k ). Thus 



Nf = Tr(U kPi 



1 



N k 



i i 



N k 



(3.88) 
(3.89) 

(3.90) 



Moreover, for each 2x2 subspace, the optimal failure probability between the two pure states 
\a k ) and \b k ) with a priori probabilities r/o and 77 1 is given by 



Qk opt 
g£opt 
Qk opt 



Vi + Vo\(a k \bk}\ 2 for J^-< \ (a k \b k )\, 
V ^0 

2y/rm~i\(a k \b k }\ for \(a k \b k )\ < J^- < 



lo + rii\(ak\bk)\ 2 for 



< 



Vo \{ak\bk) 



\{a k \b k )\ V 770' 



(3.91) 
(3.92) 
(3.93) 



In fact, the total failure probability can be expressed in terms of the canonical angles as 

1 



Q 



opt 



r . 



La 



k opt 



(3.94) 
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with for alii G [l,...,r], 



Q 



t k opt 



771 + ri Cos 2 (e k ) for 



Cos(0 k ) V Hi' 



1 / Fno 




(3.95) 



Q 



|fe opt 




(3.96) 



Q 



t k opt 



Tjo + TjiCos 2 (0 jt ) for 




<cos(e k ). 



(3.97) 



There are in conclusion numerous possible expressions (in principle 3") of the optimal failure 
probability depending on the values of the canonical angles. 

3.5.3 Unambiguous State Comparison 

Let us consider a set of n mixed quantum states {a,-} which occur with a priori probabilities {pi}. 
We are given m states out of that set and want to know with certainty whether all the m states are 
identical or not. We name this task Unambiguous State Comparison 'm out of n\ following the 
terminology introduced by Kleinmann et al. in ll28l . 

Such an unambiguous state comparison is a special case of unambiguous state discrimina- 
tion. Indeed to decide with no errors whether the m states are all identical or not, we have to 
unambiguously discriminate a first mixture of only identical states from a second mixture of non 
identical states. More precisely, we have to unambiguously discriminate the two density matrices 



where r/o = Y%=\ p" 1 and 77 1 = 1 — r/o are introduced for normalization purpose. 

In the next subsections, we are going to detail the unambiguous comparison of two pure states 
('two out of two') and a special case of unambiguous comparison of n pure states ('n out of rc'). 
We will see that those cases are reducible to some pure states scenarios and then analytically 
solvable. 

Unambiguous Comparison of two pure states 

The first case we study is the simplest situation of Unambiguous State Comparison. It involves 
only two pure states |Y+) and with a priori probabilities p + and p . We know it is 




(3.98) 



and 




(3.99) 
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always possible to write two pure states in some suitable orthonormal basis {|0), 1 1)} as |Y±) = 
05 1 0) ± /3 1 1) where 05 and /3 are real and such that 05 2 + /3 2 = 1. We can therefore denote by 
the (real) overlap between |Y+) and |Y_) as = (Y+|Y_) = 2oc 2 — 1. First of all, we write the 
two density matrices to unambiguously discriminate. Thanks to Eqn.(3.98) and Eqn.(3.99), we 
can explicitly express them as 

Po = — (pi|Y+Y + )(Y + Y + |+p 2 |Y_Y_)(Y_Y_|), (3.100) 
Vo 

pi = ^(|Y + Y_)(Y + Y_| + |Y_Y + )(Y_Y+|). (3.101) 

with r/o = p 2 + + p 2 - and r\\ = 2p + p_ so that 770 > f]i since (p + - p-) 2 > 0. Note that |YO) 
stands for |Y) <g> \ <t>). We will now show that these two mixed states are block diagonal. 

In chapter 2, we have seen that their is a freedom on the state ensemble of a density matrix. 
More precisely, a mixed state is left unchanged under a unitary mixing of its state ensemble. 
Next we remark that the density matrix pi is left unchanged if one swaps |Y+Y_) and |Y_Y+). 
Therefore, it seems natural to use a Discrete Fourier Transform to diagonalize pi. That is why, 
we can consider for pi the two unnormalized vectors 



that is to say 



l§B(: ',)(S-tP 



\b±) = X - (|Y+Y_) ± |Y_Y+» . (3.103) 



This yields the new state ensemble {—, — 1 — =, \b±)\ where 

3 V2(l±© 2 ) 



\b±) = —=^^=(\b+b-)±\b-b + )). (3.104) 
' 1 v/2(l±0 2 ) U + ' ' +/; 

We finally end up with 

pi = + © 2 )\b+)(b+\ + (1 -® 2 )\b-)(b-\). (3.105) 

It is worth noticing that, since (b + \b-) = 0, the state vectors \b±) are the eigenvectors of pi 
with eigenvalues b± = | ( 1 ± 2 ) . 
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In that form, it appears obvious that po and p\ are block-diagonal. To convince ourself, we 
simply write the different overlaps involved here. 

(Y + Y + |£+> = ; 20 , (3.106) 

2® 

(Y_Y_|Z? + > = , (3.107) 

(Y + Y + |6_) = 0, (3.108) 
(Y_Y_|£_) = 0. (3.109) 

It remains to give the optimal failure probability to unambiguously discriminate po and pi or 
equivalently the failure probability to unambiguously compare two pure states |Y±). 

In fact \b-) is orthogonal to po and to \b+) or in other words \b~) G 5^ Px fl J£p . Thanks to 
Theorem 1 1, we know that this direction \b~) can be perfectly discriminated. This direction does 
not contribute to the failure probability for unambiguously comparing |Y + ) and |Y_). We are 
left with the three dimensional subspace spanned by po and \b + ). Since po is two dimensional, 
Theorem 1 1 can again be used. It tells us that we can reduce this USD problem further and only 
consider the problem of two pure states \b+) and \b'+) with proper a priori probabilities. 

We introduce here the projection \b'+) of \b+) onto the support of po. The corresponding unit 

vector is \b+) = -^±£- cited above. We proceed as we did for the case of state filtering where here 

Ilyt is the projector onto the two dimensional subspace spanned by |& + ) and \b+). 
Theorem 1 1 tells us that the optimal failure probability Q opt is given by 



Q opt =NQ°P t (\b + },\b'[}), (3.110) 



with 



p' = ^-U^ Po U^, ^ = », iVo = Tr(pon.^) (3.111) 

p( = i-n^p 1 n^, T]( = ^, N x = Tv{ Px U^) (3.112) 

/y = /V T?o+Mm, (3.113) 

M" = {\b+},\b'[}}. (3.114) 

Let us calculate the relevant quantities N\, No and (b+\b + ). Since \b+) is an eigenvector of 
Pi, N\ simply is its eigenvalue. Thus 

1 + 2 

Ni = —- — . (3.115) 



3.5. Applications of the reduction theorems 
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To find Nq and (b'+\b+) we first have to calculate \b'L) and \b'+). We can express \b'L) in the 
non-orthogonal basis {|Y + Y + ), |Y_Y_)} of y pa so that 

\bj) = (Y + Y + |£ + )|Y + Y + > (3-116) 
20 

"Y + Y+) + |Y_Y_)). 



(l+0 2 )v/2(l+02) 
The norm of this vector therefor is 



which yields 



fr'|) = __L= =(|Y + Y + ) + |Y_Y_)). (3.118) 



+/ ^2(1 + 02) 
Since A^o = Tt(TLj^i PqII^i) = {b'\_ |po|&+) we simply obtain 



1 +0 2 

^Vo = ^— =M- (3.119) 
Finally, the last relevant quantity simply is 

~ 20 

W*+> = 11^11 = 7^02- ( 3 - 12 °) 

Considering the three possible regimes of the optimal failure probability for two pure states, 
we end up with Q° pt , the failure probability of unambiguously comparing the two pure states 
|Y±), expressed as 



20 2 fru 1 4- 2 

2(S) 2 1 4- £) 2 /rTT 90 2 

A* - *T^ + ^ fa VI^- <3 ' 123) 
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Let us note here, that we could derive the last expressions of <2 opt using the result derived for 
State Filtering (Eqn.(3.78) to (3.85)) with the following correspondences 

1 +0 2 

Po = Th—j—' ( 3 - 124 ) 

Pi = P±, (3.125) 
2@ 2 

s = ^TW (3 - 126) 

m\\ = ii^ii=t^- ( 3 - 127 ) 

In the next application of our reduction theorems to state comparison, we will use more 
properties of the Discrete Fourier Transform. 

Unambiguous Comparison of n pure states with a simple symmetry 

We propose to study the problem of comparing n linearly independent pure states | Y,-) with equal 
a priori probabilities p\ = - and equal real overlaps = (Y, | Y 7 ), Vz, j = 1, . . . ,n. 

Related to this comparison task, Eqn.(3.98) and (3.99) tell us that there is a USD problem 
that involves two density matrices po and pi and their a priori probabilities expressed as 

Po = -£|Y i ...Y ; )(Y,...Y,| (3.128) 
lo = ~~[ (3-129) 



and 



Pi = U '\ 1 A m rPo (3.130) 

n n - l -l n n - l -l 

n n - l -l 

ll = — (3-131) 
n" 1 

where 

<§ = -£P*mi- (3-132) 

Note that ^ is not a projector since the vectors |Y ; ) are in general not orthogonal. We will 
now show that these two density matrices are block diagonal and that their unambiguous 
discrimination can be reduced to n two pure states USD problems only. 

Actually we can consider the cyclic permutation C that maps |Y,) to |Y, + i) for i = 0,n — 1 
and |Y„) to |Yq) and the Discrete Fourier Transform. From now on, all the indexes are given 
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modulo n to simplify the notations. In fact, it is pretty clear that both po and pi are invariant 
under the cyclic permutation C m . We can therefore, as we have already done for the comparison 
of two pure states, use the Discrete Fourier Transform to change the state ensemble of po and 
pi. If we do so, we will see that both po and pi are block diagonal where each block is an 
eigenspace of C® n . The main reason for that is that the permutation operator C is diagonalized 
by the Discrete Fourier Transform. Importantly, the n vectors states of po are n eigenvectors of 
C®" with distinct eigenvalues (i.e. the n roots of unity). Therefore, the n vectors states of po are 
in different eigenspaces of C®". As a matter of fact, po and p\ are block diagonal where only 
one vector state of po is in each eigenspace of C m . Thanks to theorem 11 and 12, the USD of 
Po and pi is reducible to n two pure states cases. 

Now that the flow of the argumentation is clear, let us first that po and pi are invariant under 

First, we examine the action of C®' ! on Pq. 

c ®n pQC f®» = C ® B -f'(|Yi)(Y/|)®"C + ®" (3.133) 

= If c®"(|Y ; >(Y,|)®"C t ® M (3.134) 

= -V (CIY^XY^)®* (3.135) 



-£(|Y i+1 ><Y/ + i|)® B (3-136) 



n . 

i=\ 



= l£(|Y,,)(YH)®" (3.137) 

= Po (3.138) 

where the index n + 1 equals 1 since the indexes are given modulo n. We can also investigate the 
action of C the operator b, = ±£? =1 |Y,) (Y,|. 

C^C 4 = C^J^lY^YiMd (3.139) 

= -f;c|Y/)(Y/|d (3.140) 

= -£|Y^)(YH (3.141) 



(3.142) 
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Since % is invariant under C, £® n is invariant under C® n . p\ = J?\_ x £>® n — -^tzj Po where both 
£®n an( j p Q i nvar j ari t under C®", the immediate consequence is that p\ is invariant under C®" 
too. 

The Discrete Fourier Transform is the main tool of the next calculations. The matrix elements 
of U are given by 



1 



U ]k = —e^-^^, k=\,...,n. 
\ n 



The eigenvalues of C simply are the n roots of unity which can be expressed as 

Xj = e- 2i7tk ^, k=l,...,n. 
Let us briefly derive this result. In a tensor representation, C qk = 8( q+ i)k therefore 



[UCW 



'/'./ 



"Lu M c qk ul. 

qk 



1 T„(/>-l)(?-l) 

Jim- — — 



qk 



g(j-i) 



n 



2iKq 



(P-J) 



where we used the relation 



-L 



J PJ- 



(3.143) 
(3.144) 

(3.145) 
(3.146) 
(3.147) 
(3.148) 
(3.149) 
(3.150) 

(3.151) 



The unitary freedom in the ensemble of a density matrix allows us to write any density matrix 
P = JUlk\lk){pi\ as L/V/|v/)(v/| where 



v*|Vj) = £c/ yM7lMy)- 



(3.152) 



We now change the set of state ensemble of both po and p\ . In the former case, we use the 
Discrete Fourier Transform U, a (nxn) matrix acting on n non normalized vectors | Yy . . . Y/) . 

In the later case, we use the unitary transformation U on n non normalized vectors 4^|Yy) to 
change the state ensemble of ^ and therefore to change the state ensemble of p\ too. 



3.5. Applications of the reduction theorems 
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We begin with the state ensemble of po and its new a priori probabilities v,- thanks to 
Eqn.(1.6). 

y i = -I f ( 1¥ k---Y k \U* k U ij \Y j ...Yj) (3.153) 

kj 

= -L U *Uij(¥k...V k \Y j ...Y j ) (3.154) 
n kj 

= l(L u ik E u H^k • • • Y*^ . . . Yy) +£l£lfc(Y* • ■ . Y*|Y, . . . Y,-» (3.155) 

= -(®"EtftEtfy + J>fc| 2 ). (3-156) 

At that point of the calculation, two cases must be considered. On one hand there is the case 
where i = 1 and on the other hand, i ^ 1. Two properties of the Discrete Fourier Transform are 
important here. They can be summarized as 



£ \Uij\ 2 = 1 Vi. (3.158) 

The above calculation of the new a priori probabilities V; for i = 1 then leads to 

Vl = i(e"JXJ> y + 2> lt | 2 ) (3.159) 

= -(0"£^u-(v^-t/u-) + l) (3-160) 

= -(0"(VSJX-£|tf u | 2 ) + i) (3J61) 

= -(@ n (n-l) + l). (3.162) 



A similar calculation for z 7^ 1 gives 

1 



Vf = - 



(©"EOftEtfy + ElEfcl 2 ) (3-163) 

-(0 n E^(O-^) + l) (3-164) 

n * 

-(-® n Y\U ik \ 2 + l) (3.165) 
n k 

-(-0" + l). (3.166) 
n 
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Finally, po takes the form 

PD= 1 + (W " 1)0 " l<fi)<<E»il + ^-^El^XO*! (3-167) 



n it 



(3.168) 



with 



|<S>i) = — == 1 V IY/...Y,-), (3.169) 

l**> = ^ipX;^^^ |Yj • • . Y/> for / ^ 1. (3.170) 

The fundamental property of those states vector i— 1, . . . ,n is that they are eigenvectors 
of C® n with n distinct eigenvalues. Note here that C® n has the same eigenvalues than C because 
this eigenvalues are roots of unity. In other words, 

C® n \®j)=Xj\<Pj), (3.171) 

. k— 1 

with Xj = e K ~ , k= 1, . . . ,n. Indeed the operator C m acts on the vector |<E>^) as 



C m \^ k ) = C rm - 1 Ye 2in i!Lj ^ Ul \Y i ... y ¥ j ) (3.172) 
Vi-©« ■ 



./ 



£g 2i7r » - C® n \Yj...Yj) (3.173) 



Vi 



e 2,7r — » C|Y y -)®---<g)C|Y y -) (3.174) 



i y^fe^za 

Vl + (n-l)0"7 7 7 

.^f"* l e- 2l *-\Y j+1 . . . Yy+i) (3.176) 



l + (n-l)0 n y 



(3.177) 



Vl + (n-l)© n 

= g-^^IOfc) (3.178) 
= (3.179) 

By definition, po can be written in a block diagonal form where each block is an eigenspace of 

n— 1 

We follow the same technique to change the state ensemble of p\. Since p\ = n "-i_ l ^® n — 
n n-\_ l Po> we focus our interest on the matrix £, . We use the Discrete Fourier Transform U acting 
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on the n unnormalized vectors A^l^j) to change the state ensemble of £ and, as a consequence, 
of pi. 

We calculate the new a priori probabilities Vj of the new state ensemble of t, . 

Vi = -^kk\U* k Uij\Wj) (3.180) 
n kj 

= l zL U tkUij{^k\^j) (3-181) 

= l ~(L U * I fyflTO + E^tt(Y* • • ■ . . .Ty» (3.182) 

= -(©I^E^ + LI^I 2 ). (3-183) 

n k j^k k 

This calculation is similar to po's case. Only the quantity 0" is changed to 0. Therefore, we end 
up with 

v={ Kl + (»-l)®).i = l ( 3i84) 

' I J(i-0), v/^i. (3 - 184) 



Finally, ^ takes the form 



i + (»-i)e i-e 



W it 



with 



The immediate consequence is that 



Vl + (W-1)© ; 

|Y,) = -^^J^^MMlTyJfori^l. (3.187) 



' ^("-De-^^^i^^, (3 ,89) 

11 — — 



n n 1 — 1 n n 
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3. A standard form 



Moreover, the state vectors |Y ; -) of E, are eigenvectors of C therefore the state vectors |Y,i . . . Y,-„) 
of % m arc eigenvectors of C®' ! . A short calculation can verify this claim. 

CIO,) = Cj^Uflftk) 

j 

k 

= "£U jk \Y k+l ) 

k 

- L,e " \^k+i) 

k 

k 

k 



This implies that 





(3.190) 




(3.191) 




(3.192) 




(3.193) 


|TH-1> 


(3.194) 


|Y*+i> 


(3.195) 


~ l^jf ) 


(3.196) 




(j.iy /j 




(3.198) 




(3.199) 


■®A ( n|O in ) 


(3.200) 




(3.201) 



k' 

e - Mn ^\<$>j) 



C®---®C\3>ii...3? in ) = C|Oji)€ 

= A,-! | On) 



Since the state vectors of t,® n are eigenvectors of C 8 " 1 , like po, is block diagonal, where 
each block in an eigenspace of C m . The immediate consequence is that pi, linear combination 
of £® n and po is block diagonal, too. 

Let us denote Sk, the eigenspace associated with the eigenvalues A* of C®" and LT^ the or- 
thogonal projector onto Sk. We have 

po = J^pon,, (3.202) 

k 

Pi = £n,pin,. (3.203) 

k 

(3.204) 

Therefore, Theorem 12 tells us to focus our attention onto the n reduced problem defined by the 
two density matrices Pq = ^m^jk and p\ = ■ Moreover the reduced density matrix Pq 



3.5. Applications of the reduction theorems 
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simply is a pure state 

n kPo n k = \<i> k ){<i) k \ (3.205) 

since the n state vectors of po are eigenvectors of C®' 1 with distinct eigenvalues. By means of 
Theorem 11, we can reduce the USD problem of unambiguously discriminating p^ and p\ to 
the one of two pure states only. 

Finally the unambiguous discrimination of po and pi or, equivalently, the unambiguous com- 
parison of n linearly independent pure states |Y ; ) with equal a priori probabilities p\ = - and 
equal real overlaps = (Y,|Y ; ), = 1, ... ,n is reducible to n two pure states cases. 

The goal of this section was to show that the unambiguous comparison of n pure states with 
equal a priori probabilities and equal and real overlaps is reducible to some pure state case. As 
we have already indicated in the introduction, the question to know whether any unambiguous 
comparison of pure states is always reducible to some pure state cases remains opened. However, 
as expected, the unambiguous comparison of mixed states is generally not reducible to some pure 
states case 11281 . 

This concludes this chapter. In the next chapter, we will derive the first class of exact solutions 
for a generic USD problem. 
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Chapter 4 

First class of exact solutions 



The structure of this chapter is the following. In the section |4~T1 we derive three lower bounds on 
the failure probability to unambiguously discriminate two density matrices in three regimes of the 
ratio between the two a priori probabilities. Our derivation uses the Cauchy-Schwarz inequality 
and allows us to look for necessary and sufficient conditions to reach the lower bound in each 
regime of the a priori probabilities. In section l4~2l we report the notion of parallel addition that 
leads to some useful relations for USD in connection with our first reduction theorem. In section 
4.3, we finally derive the main result of this chapter as a theorem: a necessary and sufficient 
set of two conditions for the failure probability to reach the bounds are given. We also give the 
corresponding optimal POVM. 

With that result, we give the optimal USD POVM of a wide class of pairs of mixed states. 
This class corresponds to pairs of mixed states for which the lower bounds (one for each of the 
three regimes depending on the ratio between the a priori probabilities) on the failure probability 
Q are saturated. This class in nonempty since it contains some pairs of generic mixed states 
as well as any pair of pure states. For those pairs, we provide the first analytical solutions for 
unambiguous discrimination of generic mixed states. This goes beyond known results which are 
all reducible to some pure state case as we have seen in chapter 2 and 3. 



4.1 Lower bounds on the failure probability 

The failure probability Q of a USD strategy is given by Q = Qj, where Qj = r\[Yx{E7pi) . From 
this definition we immediately see that Qj < r]j. In this chapter, we consider the USD of two 
signal states po and pi that are mixed states with a priori probabilities T]o and r\\. Accordingly, 
our POVM contains three elements {Eq,E\,E?} which correspond respectively to the conclusive 
detection of po, to the conclusive detection of pi and to an inconclusive result. The failure 
probability then equals Q = Qo + Q\. 
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4. First class of exact solutions 



Our interest is first focused on the product QoQ\. We can give a lower bound expressed in 
terms of the fidelity F = Tr ( «J sfpop 1 y/po) of the two states po and p\. The bounds, formulated 
in the following theorem, are tighter than those given in chapter 2. Moreover, we pay additional 
attention to the condition under which the bounds can be reached. 

Theorem 13 Lower bound on the product QqQ\ 

Let po and p\ be two density matrices with a priori probabilities r\$ and r\\. We define 
the fidelity of the two states po and pi as F = Tr( \J y/pop\ \7po)- Then, for any USD 
measurement with inconclusive outcome E?, the product of the two probabilities Qq and Q\ 
to fail to identify respectively the state po and pi is such that 

QoQi>lomF 2 . (4.1) 
The equality holds if and only if the unitary operator V arising from a polar decomposition 

VPoVPi = \JVPoPi\/PoV (4.2) 

satisfies 

V t v /poV^= ayfplx/Wi (4.3) 

for some a G 1R + . 

Before we turn to the proof of this theorem, note that relation (14.31) implies a condition on 
the optimality of a USD POVM J32l[33l|45l. It is clear that optimality of a specific USD mea- 
surement implies that the conditional states after the inconclusive results do not allow further 
USD measurements as we already discussed it in chapter 3. This condition is satisfied, for exam- 
ple, when the supports of the conditional states coincide. We find a stronger property whenever 
equality holds in Theorem 13. Indeed, if we have V* y/po^/E? = a^/pi^/E? with a G 1R + , then 
it follows immediately that y/Ejpoy/E? = a 2 y/E?piy/E?. This means that the conditional states 
corresponding to inconclusive results must be identical up to normalization. Therefore no infor- 
mation whatsoever about the signal state can be extracted from these conditional states. 

Proof of Theorem 13 This theorem was stimulated by the proof of the nonbroadcasting theo- 
rem [|46l . The basic ingredient for the derivation of the bound is the Cauchy-Schwarz inequality: 

Theorem 14 WTi Cauchy-Schwarz inequality 

Ifx andy are members of a unitary space then \\x\\ \\y\\ > \(x,y)\. 

The equality holds if and only ifx = ay for some a in C. 

A unitary space is a complex linear space y together with an inner product from J? x J? to C. 
Therefore the complex space of bounded operators acting on a Hilbert space is a complete unitary 



4.1. Lower bounds on the failure probability 
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space (i.e. every Cauchy sequence converge) if we consider for two elements A and B the inner 
product Tr(A5^). The Cauchy-Schwarz inequality then takes the form ^Tx{AA^)yjTx{BB^) > 
|Tr(A5 t )| where equality holds for A = aB, a in C. 

Let us now consider a POVM element E k and two density matrices po and p\. We can 
decompose these three operators as pi = \/pT\/pi and E k = \[E\-jE{ and po = v /pot/f/^y / po 
where U is an arbitrary unitary transformation coming from the freedom in the decomposition 
of a positive semi-definite operator. Hence we obtain from the Cauchy-Schwarz inequality with 
A = U^^^E~ k and B = Jp[jE k 

^Tr(E kPo )^Tr(E kPl )> |Tr(^ ^pV^WpT)! = (Tr^^P^/pI)!- ( 4 - 4 ) 

By Theorem 14, the equality holds if and only if U* \f P Q^/E k = a^fpl^/El, for some a E C. 
We now consider a USD POVM {£jt}jt=o,i,?. Using the fact that Tr(E Pi) = Tr(£ip ) = 0, we 
find for Eq and E\ 

0=^Tr{E oPo )^Tr{E oPl )> |Tr(t/ t v /p^ 0v / PT)l. ( 4 - 5 ) 
0= ^Tr{E lPo )^Tr(E lPl ) > \Tr(U^E 1 ^pT) |. (4.6) 

This simply means that Tr(U' JPqEq-s/Pi) = Tr{U^ ^/ P qEi^J P ~\) = 0. For E?, we obtain 

^Tr(E 7Po )^Tr(E ?Pl )> \Tx{U^ ^E ly fpl)\. (4.7) 

From this it follows that we can write 

^Tr{E 7Po )^Tr{E ?Pl )> \Tr(U*JtoE ? Jfr) + + 0| = |Tr(tfVPoVPi)l » ( 4 - 8 ) 

where we used the relation Y*k.Ek = 1. Furthermore, the inequality (14.81) must hold for any unitary 
matrix U so that we find 

y/ Tr(£ ? p ) yj Tt{E ?Pi ) > max \Tr(U* y/pSy/pt) | . (4.9) 

Here, again, the equality holds if and only if a unitary operator U max which maximizes the right 
hand side satisfies 

^maxV / P^V / ^= <xVP~lVE~? (4-10) 

for some a E C. To find the unitary matrices U max that maximize |Tr(£/ ' v^Ak/Pi) | we use the 
following lemma: 
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4. First class of exact solutions 



Lemma 2 For any operator A in the space M„ ofnxn matrices we find 

max\Tr(AW)\ = Tr(\A\) (4.11) 
w 

where the maximum is taken over all unitary matrices. The maximum is reached for any unitary 
operator W that can be written as W = V^e l< ^. Here e 1 ^ is an arbitrary phase while the unitary 
operator V is defined via a polar decomposition 

A=\A\V (4.12) 

with \A\ = Va~A^ = VVA^AV^. (See proof in Appendix B.) 

Let us introduce the operators Fq := | v /poy / pTl — s/y/popiy/po and F\ = V^FqV = 
■J v /pTpoy / pT, which are motivated by the polar decomposition 

sfpoVpi=F Q V = VF l . (4.13) 

These operators are related to the fidelity of the two density matrices through the relation F = 
Triy/y/pipoy/pi) = Tr(F ) = Tr(Fi) ED. Thanks to lemmad Eqn. implies 

^Tt(E iPo )^Tt(E iPi ) > |Tr( I v^v^TI ) | = Tr(| y/fry/pi\ ) (4.14) 

where equality now holds if and only if U max in (14.101) arises from a polar decomposition of 
^/po^/p~l■ In other words, we have 

y t e^ v ^ v ^= a^fpl^Ej (4.15) 

for some aeC. 

Next we use the definitions of the partial failure probabilities <2, = r\{Ix{E?pi) and 
choose the phase e'^ to be the same as the phase of a in (14.151) to obtain the desired 
inequality <2o<2i > loliF 2 . Equality in the previous equation then holds if and only if 
y^y/PoVE? = a^fply/Ei, for some a e R + . This completes the proof. ■ 

We can now derive the bounds in the different regimes of the ratio ^- between the two a 
priori probabilities. Actually, the procedure is to find the minimum of the failure probability Q = 
Qo + Qi under the constraints of the previous derived inequality <2o£h > loliF 2 . According to 
Theorem 13, we can provide the necessary and sufficient condition for equality. 



4.1. Lower bounds on the failure probability 
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Theorem 15 Lower bounds on the failure probability 

Let po and p\ be two density matrices with a priori probabilities r/o and r\\. We define 
the fidelity F of the two states po and pi as Tr{*J ^/p~oPi y/Po) ■ We denote by Pq and P\, the 
projectors onto the support ofpo and pi. Then, for any USD measurement with inconclusive 
outcome E?, the failure probability Q obeys 

Tr{Pip ) 



Q>2^m~iF for ' <JJ±< (4.17) 




fi>*to T /p n s +VMPoPi) for . <Wgi. (4.18) 
Tr(Popi) Tr(Popi) \ T] 

Equality holds if and only if the unitary operator V arising from a polar decomposition 
y/Poy/pi = y^/PoPu/Poy satisfies y/poy/E? = a^fpl^fEi, with a = Tr ^ Po \ a = 
and a = Tr (p oPl ) in the the first, second and third regime, respectively. 



Proof First of all, according to Theorem 13, we know that for any USD measurement the 

Qo 



inequality QoQ\ > t}qT]\F 2 i.e. Q\ > riol ^ lF holds. It follows that the failure probability is lower 



bounded as 

e>2o+^. (4.19) 
Qo 

Since we are interested in a lower bound on Q, let us consider the case where equality holds 
in Eqn. (14.191) . In this case, we have 

GoGi = rio^F 2 (4.20) 

Q = Q + ^fl (4.21) 
Qo 

From Theorem 13 we know that Eqn. (14.20l) holds if and only if '-y/poVE? = a^fpl^fEj, for 
some a G R + . We will now connect a to the other quantities. The previous relation implies, via 
the respective definitions, that 

Qo = a 2 -^£i . (4.22) 

r?i 

The former relationship corresponds to the proportionality between two vectors of the vector 
space of bounded operators while the latter relationship corresponds to the proportionality be- 
tween their norms. We can combine the two equations (14.201) and (14.221) to 

Qo = arioF (4.23) 

Qi = ^-llF. (4.24) 
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4. First class of exact solutions 



So the final statement is that Q = Qo+ r,0 ^ ) lF if and only if ^/poy/E? = a^/piy/E?, where a 



Qo 

now is explicitly related to the other parameters as Qq = ar\$F and Q\ = ^f]\F '. 

Second, we have to derive a range constraint on Qq and Q\. We know already that Q\ < 
r\i. Moreover, from work by Herzog and Bergou in |29], we learn that r\oYx(P\po) < Qq and 
7]iTr(i 5 oPi) < Q\ . Indeed, from the structure of the USD POVM elements, we have E + E\ + 
£7 = 1 with ,5^Eq C Jf Pl and 5?E\ C J£p - We consider only the non-trivial case where the 
supports of po and p\ are not identical. Then the structure must be such that E\+E? = P\+R 
where P\ is the projection onto the support of p\ and R is a positive semi-definite operator with 
support S^r C Jff Pl which satisfies Eq + R = P^- otherwise Tr(.EoPi ) 7^ 0. Then it follows that the 
partial success probability Pq is Pq = 7]oTr(£'oPo) = T ?oTr(i , j 1 po) — 7]oTr(7?po) . In our non-trivial 
case we will have Tr(i?po) > as soon as R 7^ 0. This yields Pq < r]oTr(P^ po) or equivalently 
Qo > ^oTr(PiPo). In the same way, one can find Q\ > T]iTr(Popi)- We then have 

TfoTr(PiPo) < Qo < t?o, (4.25) 
T?iTr(^pi)<Qi<Tji. (4.26) 

These two constraints can be combined in 

F 2 

77oTr( J P 1 p ) < Qo < ?7o T (T3 - v (4.27) 

This can be seen as follows. Since Q\ = 7,0 , the constraint (14.261) on Q\ takes the form 

F 2 

rioF 2 <Qo< Wo— r-. (4.28) 

We now have two lower bounds and two upper bounds on Qq ( (14.251) and (14.281) ) and we want 
to find the tighter ones. To do that, let us consider the USD POVM given by {E*? — P\,Eq — 
P^,E\ = 0}. Thank to Theorem 13, we find f]oTiiF 2 < T7oT]iTr(Pipo)Tr(Pipi) or in other words 
r\oF 2 < T] Tr(P 1 p ). We can also consider the USD POVM given by {E? = P ,E = 0,Ei = 
P^} and with Theorem 13, we have t/o Tr ^ )Pi ^ < Wo- Finally, we obtain rioTr(Pipo) < Qo < 

Next, we define the function q (Qo) = Qo + 770 q iF and minimize it under the constraint 



Qo 

T]oTr(/ , ipo) < Qo < r/o Tl .^p QPl s ■ The resulting minimum will constitute a lower bound for Q. 
The function q(Qo) is convex (4S(Qo) > 0) and, therefore, it takes its minimum at the point 

"Qq 

2™ n where the derivative vanishes (Q ) = yielding Q^ m = ^/r]oV^F) or at the limits of the 
constraint interval (<2™ n = r}oTr(P\po) and <2™ n = lO j^p^ )- That gives us the minimum of the 
function q(Qo) in three different regimes. In the first regime we have q m m(Qo) = T?oTr(Pipo) + 
7]! and 2™ n = T 7o Tr( J P 1 p ) if ^WTiF < T 7o Tr( J P 1 p ) that is to say if ^ < 
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f < fnL 

Hr(fl)pi) - V rj - 

As a result we obtain lower bounds for the failure probability 2 in three regimes as given 
in Eqn. (14.161) . For each regime, the value of Qq which minimized q(Qo) is given and via 
Eqn. (14.231) we find the corresponding value that a has to take. We read off the values as 
a = Tl ( p iP°\ a = and a = Tr ^ pi ) f° r the first, second and third regime, respectively. This 
completes the proof. ■ 

Let us note that, by construction, those bounds are tighter than the ones in chapter 2 ll26l . 
Indeed, one could recover the three bounds in ll26l by looking for the minimum of the function 
q(Qo) under the weaker constraints rjoF 2 < Qq < r/o as we will show in the last section of this 
chapter. 



Before deriving the central theorem of this chapter and then provide the first class of exact so- 
lution for USD of two generic mixed states, we will first recall some useful results of linear 
algebra. We denote by the pseudo-inverse of a matrix M, which has not necessarily full 
rank. The pseudo-inverse can be defined via the singular- value decomposition of M = UDV as 
= UD l V, where U and V are unitaries and D is a positive semi-definite and diagonal ma- 
trix. Whenever M is of full rank, the pseudo-inverse coincides with the inverse. In general, it is 
not known how to express the pseudo inverse of a sum (A + in terms of the pseudo inverses 
A -1 and B 1-1 Il48ll4~9~ll . However, a related operation A (A + B) _1 B, called parallel addition and 
denoted by A : B has been defined and studied in 1969 by Anderson and Duffin and will turn out 
useful in our context. 

Lemma 3 ll48l Let A and B be two positive semi-definite matrices in M n , then the support S^a-.b 
of A : B is given in terms of the supports of A and B as 



(See proof in Appendix C.) 

Next let us recall the first reduction theorem for USD of mixed states (Theorem 9). We 
consider the problem of discriminating unambiguously two density matrices po and pi with a 
priori probabilities 770 and T]i . We denote by ro the rank of po and by r\ the rank of pi . A general 
USD problem can satisfy ro + r\ > d, where d is the dimension of the Hilbert space J$? spanned 
by the two states. This means in particular that the two supports can overlap. 



4.2 Parallel addition p E 1 



Pi 



yA:B = y A ny B . 



(4.29) 
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4. First class of exact solutions 



In the first reduction theorem it has been shown that any such USD problem can always 
be reduced to the one of discriminating p' Q and p[, two density matrices of rank r' Q and r\ 
with a priori probabilities T]q and rj[, spanning the same Hilbert space of dimension 
d = r' + /j . Indeed we can split off any common subspace of the supports y po fl y pi to end 
up with y p i fl y p i = {0}. As we have already seen, two supports do not overlap if and only 
if rank(pQ) + rank(p[) = rank(p' + p[) holds. In such a reduced case, Lemma [3] implies 
y p ': P ' x = that is to say p' Q :p[= 0. 

We defining E := p' Q + p[ to write the parallel addition as p'^L~ l p[. Since rank(p' Q + p[) = 
dim(jtff), we end up with E having full rank and EE -1 = 1j%> . We therefore have the following 
corollary to Lemma [3l 

Corollary 3 Let po and pi be two density matrices spanning a Hilbert space M'. Let E be the 
full rank operator defined as the sum of these two density matrices. 

If y p , ny pi = {0} then poE-Vi = 0. 

According to the first reduction theorem we can, without loss of generality, consider only 
USD problems of two density matrices without overlap of their supports. In the following, we 
consider two density matrices po and p\ (which are positive semi-definite matrices) such that 
y pi) fl y pi = {0} or equivalently rank(po + p\) = rank(po) + rank(pi) = dim(j^). As ex- 
plained above, for such a problem, poE _1 pi = 0, with E = po + pi having full rank. This leads 
to 

Pi = PiTr l p h i = 0,1 (4.30) 

since EE -1 =\^. The projectors onto the supports of p 7 , i = 0, 1, can then be written as 

Pi = v^E- 1 ^, i = 0,l (4.31) 

To finish, let us precise that the two density matrices involved in a standard USD problem 
fulfill all the above properties since they do not overlap. 

4.3 Necessary and sufficient conditions - first class of 
exact solutions 

We are now ready to derive the main result of this chapter. The first part of this result gives 
compact necessary and sufficient conditions for a pair of mixed states to saturate the bounds of 
the failure probability Q. The second part gives the corresponding POVMs in an explicit form. 



4.3. Necessary and sufficient conditions - first class of exact solutions 
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Theorem 16 Necessary and sufficient conditions to saturate the bounds on the failure probabil- 
ity 

Consider a USD problem defined by the two density matrices po and pi and their respective 
a priori probabilities r\o and r\\ such that their supports satisfy y pQ fl S^ Pl = {0} (Any USD 
problem of two density matrices can be reduced to such a form according to Theorem 9). 
Let Fq and F\ be the two operators ^/^/popu/po and v /pTpo^/PT- The fidelity F of the 
two states po and pi is then given by F = Tr(Fo) = Tr{F\). We denote by Pq and P\, the 
projectors onto the support of po and p\. The optimal failure probability Q opt for USD then 
satisfies 



Q 



opt 



m 



+ T]oTr{P { po) & 



Po _mpo) Fo > 



Q 



opt 



Q 



opt 




pi 



Tr{PiPo 



■Fx >0 



J ^ F °-° f 0r T < P ^ < 



;F >0 



^Fi>0 



for 



1 < Tr(P lPo ) 





for 



Tr(P Pi 
F 



Tr(P Q p h 



V t?o 
(4.32) 



The POVM elements that realize these optimal failure probabilities, if the corresponding 
conditions are fulfilled, are given by 



Eo 



po (po - aFo) 1 
Way/po~+ -^=y/p[V*^j F ^vWPo + ~^ V ^P l 



(4.33) 



with a 



T '( p i p °} far the first regime, a 



^ for the second regime and a 



TriPoPx) 

for the third regime and where the unitary operator V arises from a polar decomposition 



Proof of Theorem 16 First, we give a proof for the necessary conditions. 



Proof for the necessary conditions From Theorem [15] we know that the bounds on the failure 
probability are satisfied whenever V^^/poE? = a^/piE? with a = Tl ^ Po \ a = and a = 

j r ^P oPl ) f° r the three regimes, respectively. 
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4. First class of exact solutions 



We replace F? by 1 — Eq — E\ , multiply on the left by V and on the right by y/po- This leads 
us to 

po - aF = y/foEoy/fo (4.34) 



where we used the relation (14.131) ^/po^fpl = FqV and the fact that the support of p,- and Ej are. 
orthogonal for i ^ j. Indeed, in Lemma 1 , we have seen that Tr(F,-py ) = <^ Ejpj = because F, 
and pj are positive semi-definite operators. The right hand side in (14.341) is positive semi-definite 
because of the form AA^ with A = ^/p^/Eo. Thus po — OcFq must be positive semi-definite as 
well. A similar calculation where we multiply on the right by Jp~[ instead of by J~pQ leads us to 



pi-^Fi = ^p[E ly fp[ (4.35) 
which is again a positive semi-definite operator. 

With that we have proved that if equality holds in the bounds of Theorem [T5] then we have 

Po - uF > 



(4.36) 

Pi-5*i>0 

which form, therefore, necessary conditions for equality in the bounds of Theorem [T5l 

Proof for the sufficient conditions Now we start with the assumption that the conditions 
(14.361) are fulfilled. Then we can construct an explicit POVM saturating the bound, therefore 
providing that the conditions are sufficient. Let us define the following POVM elements : 

E = E- 1 v ^(po-aF ) v ^E- 1 (4.37) 



( pi - ) yfplYT 



First, let us verify that this is indeed a valid POVM. The three operators are positive 
semi-definite since they are of the form A^MA where M is a positive semi-definite operator. In 
the first two cases this is true because of the conditions (14.361) . in the third case it follows from 
the positivity of Fo. The three operators sum to identity, Fo + E\ + E? = 1, as can be checked by 
straightforward though lengthy calculation, making use of Eqn. (14.131) . Next, we have to check 
that the given POVM is a valid USD POVM, that is, Tr(p Fi) = Tr(piF ) = 0. This relation 
holds since the supports of po and pi do not overlap. Therefore, corollary 3 applies and we 
have poE _1 pi — from which follows that y/pTjE^pi = and y/p\Hr l pQ = 0. Finally, one 
can check in a straightforward though lengthy calculation, exploiting the properties used in the 
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previous checks that this POVM leads to the three desired failure probabilities. This completes 
the proof. ■ 

Let us first note that the assumption about the non-overlapping supports was only used to 
prove the sufficiency of the conditions. Their necessity does not require this assumption. 

Moreover given a pair of two density matrices with their a priori probabilities, the middle 
regime does not always exists. A necessary condition for its existence is 

Tr{P lPo )Tr(P pi)<F 2 (4.38) 

as pointed out by Ulrike Herzog in ll29l . 

To conclude the presentation of our first class of exact solutions, we would like to repeat that 
only the first reduction theorem is needed to derive Theorem 16. In chapter 6, we will provide 
pairs of density matrices that fall in this class as well as pairs of density matrices that are not 
included in it. It means that this class contains pairs of density matrices but does not cover all 
pairs. 



4.4 The two pure states case revisited 

It is possible to use Theorem 16 for two pure states |Y±). We change here the label of the 
two states from '0/1' to ' + /— ' since one can always write two pure states |Y±) = a\0) ± 
/3|1) where a and /3 are real and such that a 2 + /3 2 = 1 in some suitable orfhonormal basis 
{|0), 1 1)}. For two pure states, the operators F± are easy to explicit. Indeed F + = F\Y+) 
andF = F|Y_)(Y_| withF = |(Y+|Y_)| = \2a 2 -l\. Moreover one has the simple relation 
Tr(P+p_)=Tr(P_p + )= J F 2 . 

The conditions in Theorem 16 then take the following form: 



[\-F 2 )p + >0 for J—<F 



(4.39) 



for F < J— < - 



1- 




1+ 



1 



(l-^)p_>0 for-< A /^ 



n- 



Since (1 — \p^F) an d (1 — \J~^~F) f° r F — \f^Z — ^ ran § e between and F 2 , the con- 
straints above are always fulfilled and our result reduces to that of Shimony and Jaeger. Moreover 
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4. First class of exact solutions 



we can give the POVM elements in a compact form thanks to the operator E 1 . The choice of 
our basis yields 



a 2 ±a/3 
±aj8 /3 : 



P±= ( , „o o2 ) ( 4 - 4 °) 



such that 



1 / a" 2 
"2 V j8" 



£ =«( „ p-2 )• ( 4 - 41 ) 



It is therefore easy to write the optimal USD POVM as follows 

,-2 l 




a 2 ^6 



a/3 

and 



«fl P 



(4.42) 



(4.43) 



1 -£+-£_ (4.44) 
— for the second regime and a = A 



with a = F for the first regime, Of = A / -jj- for the second regime and a = A for the third regime. 



This expression of £± leads naturally to the desired failure probability Q opt = F(ar}+ + ^-) 
with the respective as. 

We can go beyond this remark and investigate under which conditions our bounds reduce to 
those given in chapter 2. The bounds derived by Rudolph et al. in [|26Tl take the form 



2 opt > iH + VoF 2 for x f^<F, (4.45) 

V 

2 opt > 20p7TFforF< A /^-<4> 

V T70 F 



2 1 / /m 




e opt > no+m^ for-<, 

^ V 

Actually one can find Rudolph's bounds following the argumentation in the proof of Theorem 
15 but using the weaker constraint rjoF 2 < Qo < T}q. This means in particular that our bounds 
are tighter. To convince ourself, we can nevertheless consider our bounds and Rudolph's bounds 
in the five regimes of the ratio given by 

< P <m&<r <l (4.46) 
F Tr(Popi) F 
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Note that this ordering is due to the Theorem 13 that tells us that F < r ^ P() ^ since F 2 < Tr(Pipo) 

and Tr (^p \ < j since < Tr(PoPi)- On the other hand, the inequality Tr (^ P() ) < frT^pTj * s 
not always fulfilled as we already discussed (see Eqn. (4.38)). We can now compare the two 
bounds in each regime. 



In the middle regime given by Tr (^ p °) < < Tr(i'oPi), the two bounds are equal. In 



the second regime given by F < \J~^ < Tr ^ Po \ Rudolph's bound still equals the overall lower 
bound 2 v /rjo7|j7 7 and is therefore less or equal than our bound. In the third regime given by 
Tr (P {)P] ) < \J~^ < j?, a similar argument holds: Rudolph's bound still equals the overall lower 
bound 2^/ffoff[F and is therefore less or equal than our bound. 



In the outer regimes, things are a bit more subtle. We must again consider the function 
<l(Qo) — Qo + ^go — ■ This function decreases for < Qo < y^f\^r\\F and increases for 
VWhF < Qo- 

In the first regime, we have by definition y^j- < F < Tr( ^| Po ) (See Eqn. (4.38)). We can multiply 
this inequality by T}qF to get ^rjoTfiF < 7]qF 2 < r\QTx(P\po). For that range, q{Qo) increases 
so that Q{t)qF 2 ) < <2(noTr(Pip )) or in other words: T) Q F 2 + t} x < ri Tr{P lPo ) + r)x^ 



PoJ 



In the fifth regime, we have Tr ^ pi ) < j < \J~^ (See Eqn. (4. 50)). We can again multiply this 
inequality by rjoF to get lo^^^j < ?7o < ^/rjoliF. For that range, q(Qo) decreases so that 
Q^t^pT) ) ^ 2(^0) or in other words: Voj^Tj + ^i Tr (A)pi) < lo + I1F 2 . 



Since our bounds are tighter, Rudolph's bounds are reached if and only if, first, the conditions 
in Theorem 16 are fulfilled and, second, the equalities Tr(PoPi) = Tr(Pipo) = F 2 hold like in 
the pure state case. Let us now state the corresponding theorem and give the only part of the 
proof that changes with respect to theorems 15 and 16. 
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4. First class of exact solutions 



Theorem 17 Necessary and sufficient conditions to saturate the bounds in K2M 
Consider a USD problem defined by the two density matrices po and pi and their respective 
a priori probabilities T]o and r\\ such that their supports satisfy y pQ fl S^ Pl = {0} (Any USD 
problem of two density matrices can be reduced to such a form according to Theorem 9). 



Let Fq and F\ be the two operators ^/popi^/po and y^/pTp^/pT- The fidelity F of the 
two states po and pi is then given by F = Tr(Fo) = Tr{F\). We denote by Pq and P\, the 
projectors onto the support of po and p\. The optimal failure probability Q° pt for USD then 
satisfies 



)0 pt_ „ r 2 _ J Po-FF >0 




The POVM elements that realize these optimal failure probabilities, if the corresponding 
conditions are fulfilled, are given by 




E = E- 1 v ^(po-aF ) v ^E- 1 (4.48) 
E X = Z-^x^Pi-^Fxj ^plZ- 1 

1 ( . fry /TTTj /TvTT/t \ El. I , frt /7T7 J _T7 /TTTIV^ 1 



Z- 1 [ VaVPo~+^VP~iV J \F \VaVPo~+^V^p- l )Z 



with a = F for the first regime, a = J ^ for the second regime and OL = j for the third 
regime and where the unitary operator V arises from a polar decomposition ^/Po y/pT = 
V 'y/p6piy/p6V : 

In the first regime a = F implies that E\ = 0. The resulting POVM has to be a projective 
measurement with projections onto the support of pi and onto its orthogonal complement, i.e. 
Eq = Pj 1 , E\ = and E? = P\. A direct proof from the explicit expressions in Eqn. (4.48) is 
difficult, however a simple reasoning allows to verify this statement. We consider only the non- 
trivial case where the supports of po and p\ are not identical. Of course, a two-element USD 
POVM satisfies Fo + F? = 1 with S^e C S? Px . Then its structure must be such that F? = P\ + R 
where Pi is the projection onto the support of pi and R is an operator with support S^r C <W Pl 
which satisfies F + R = P±. Then it follows that Q = 7]i + 77 Tr(Fip ) + 7] Tr(Fpo). In our 
non-trivial case we will have Tr(Fpo) > as soon as R ^ 0. Therefore we find as an optimal 
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solution within this class of two-element USD POVM, the POVM with R = leading to E? = 
Pi and Eq = P^~. We can actually write the failure probability as <2° pt = + VoF 2 . Indeed 
pi = jF\ then p\ = jv^/p[pQ^/p\■ This implies F 2 pi = P\poP\ and finally Tr(Pip ) = F 2 . 
This is consistent with the results derived above and gives the correct failure probability. In the 
third regime, we have a = j. Therefore Eq = and the corresponding POVM is a projective 
measurement with Eq = 0, E\ = Pfj-, E? = Pq. 

Proof of Theorem 17 We will only derive the three minima of the function q(Qo) = Qo + 
^gjj- since the remaining part of the proof does not change (the proof correspond to Theorem 
15 where the bounds are derived). Here we consider weaker range constraints on <2o an d Qi '- 
< Qo < 7]o and < Q\ < T]i. We then minimize q(Qo) under the constraint rjoF 2 < Qq < T}q. 
Again, the function q(Qo) is convex (-rz&(Qo) > 0) and, therefore, it takes its minimum at the 



point <2™ n where the derivative vanishes (j^(Qq) = yielding <2™ n = \Jt]qT]\F) or at the 
limits of the constraint interval (<2™ n = t\qF 2 and Q™ m = t]q). That gives us the minimum of 
the function q(Qo) in three different regimes. In the first regime we have ^min(Go) = VoF 2 + r\\ 



and <2™ n = ^o^ 2 if \Jr\or\iF < t\qF 2 that is to say if J^- < F. In the second regime we 



have <7 min (£ ) = 2 V / W?T^ and Q™ m = ^/fjWiF if F < * f < I. The third regime gives 



As a result we obtain lower bounds for the failure probability Q in three regimes as given in 



Eqn. (14.471) . Since Qq = ccrjoF, we read off the values of a as a = F, a = , ^- and a = j for 



In the next chapter, we will derive a second class of exact solutions. This class is concerned 
with pairs of geometrically uniform states in four dimensions. 






the first, second and third regime, respectively. This completes the proof. 
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Chapter 5 

Second class of exact solutions 



In this chapter, we derive three important results. First we derive a theorem concerned with the 
rank of an optimal USD measurement. Next, we propose a corollary which is interested in the 
spectrum of an optimal USD measurement. Finally we give the main result of this chapter, a 
second class of exact solutions. This class corresponds to any pair of geometrically uniform 
states in four dimensions. To be proved, this result requires most of the theorems previously 
derived in this thesis. 



5.1 Overall lower bound and rank of the POVM ele- 
ments 

The maximum rank r'™ x of a USD POVM element Ej, i = 0, 1 is 

r^ x = dim{j(r pi ), (5.1) 
r^ x = dim(jep ). (5.2) 

(5.3) 

In the case where S? Po PI 5? Px = {0}, the maximum rank of the USD POVM elements i = 0, 1 
is dim(y pi ), the rank of the mixed states p,-. Indeed Ei has support in J(f Pj , i,j = 0, 1, j ^ i and 
therefore, if ^ Po n y pi = {0}, 

rank(Ej) < dim(jC Pj ) (5.4) 

< dim(JF)-dim(y Pj ) (5.5) 

< dim(y pQ ) + dim(y pi ) - dim{y pj ) (5.6) 

< dim(y p .), i = 0,l. (5.7) 

Note that in Chapter 3, we already proved that 

r™* = min(dim(y po ),dim(y pi )). (5.8) 
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The first theorem of this chapter states that the two POVM elements Eq and E\ of an optimal 
USDM both have maximum, rank only if the two operators po — \J~^Fq and p\ — \J~^F\ are pos- 
itive semi-definite. The attentive reader can recognize the two operators involved in the middle 
regime of Theorem 16. 

Theorem 18 Rank o/Eq and E\ 

Consider a USD problem defined by two density matrices po and pi and their respective a 
priori probabilities tt\q and r\\ such that their supports satisfy 5^ Pq fl 5 /fi Pl = {0} (Any USD 
problem of two density matrices can be reduced to such a form according to Theorem 9). 
Consider also an optimal measurement {Eq PI ,E opt ,E opt } to that problem. Let Fq and F\ be 
the two operators y^/popu/po and ^/pipQ^/pi. The fidelity F of the two states po and 
pi is then given by F = Tr(Fo) = Tr{F\). 

If the two POVM elements Eq P and E opt have maximal rank dim(j7 ? Po ) and dim(y pi ), 
respectively, then 




(5.9) 



Proof We consider an optimal measurement for unambiguously discriminating two mixed 
states po and pi. We can therefore use the necessary and sufficient conditions derived by El- 
dar OBI . We recall them here. Necessary and sufficient conditions for a measurement {E{\, 
k =?,0, 1 to be optimal are that there exists Z > such that 

ZE? = 0, (5.10) 

£o(Z-TfoPo)£o = 0, (5.11) 

£i(z-77ipi)£i=0, (5.12) 

Pt(Z-r]oPo)Pt>0, (5.13) 

P - L (Z-i]ipi)^6 L >0. (5.14) 

If Eq and E\ have maximum rank and Eqn. (5.11) and Eqn. (5.12) are fulfilled then the two 
Hermitian operators P^(Z — rjopo)P^ and P 1 (Z — tjiP^Pq 1 - must vanish. Indeed the situation 
is the following. We consider two positive operators A and B, with A full rank and ABA^ = 0. 
We can see this relation as of the form CC^ = with C = A^/B. Moreover, such an equation 
CC^ = is equivalent to C = for any matrix C (See Appendix A for a proof of this statement). 
Consequently, ABA^ = is equivalent to A\^B = 0. Finally, since A is full rank A -1 exists and 
B must vanish. 



5.1. Overall lower bound and rank of the POVM elements 
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In Eqn. (5.11) and (5.13), we have A = E and B = P^(Z - 7] p )P^-. In Eqn. (5.12) and 
(5.14), we have A — E\ and B = P^(Z - Tjipi)P^\ As a result, P^(Z - r] p )P^ and P^{Z - 
r\\p\)PQ must vanish if Eq and E\ have maximum rank. Finally to prove the statement of the 
theorem we can show the following equivalence: 

f ,, ZE? = ° , f P0-\/fo F ^° 

3Z > such that I Pq{Z- t]\P\)Pq = ^{ V r 1 (5 ' 15) 

I p l Hz- Wo )pt = o ( Pi-V^° 

where Pq{Z — x\\P\)Pq and P^{Z — t\qPq)P^ are positive semi-definite operators. To prove 
this statement, we proceed by equivalence. 

Since the two supports do not overlap, we can make use of the full rank operator XT 1 = 
(po + Pi)~ ! introduced in chapter 4. Let us repeat here that its main property is 

p i Yr i pj = p i 8 i] , r = 0,l. (5.16) 
As a consequence, we get the interesting relations 

PoE- 1 = Po^Pt, (5.17) 
P^po^- 1 = Pt- (5.18) 

Indeed po^ 1 = p £ _1 (^i + p t) = Po^PiPf 1 + Po^ 1 ^) = po^P^- Moreover, if = 
Pj 1 ! = P^(po + pi)£ _1 = P^PqL~ 1 . The same relations are of course true when we swap 
and 1. 

PlE" 1 = piZ- 1 /^, (5.19) 
/^PiE- 1 = Pi~. (5.20) 

It follows that the two equalities P^(Z - rj p )P^ = and Pq(Z - T^pi)^ = are equivalent 
to poYr l {Z- r\Qpo)Yr l pQ = and p{L~ l {Z- r\ip\)Hr l pi = 0. Hence the assertion 



ZEj = 

3Z > such that <( P^{Z-T]ip l )P^ = Q (5.21) 
if (Z-TfoPo)if = 



can be replaced by 



7F? = 

3Z > such that <J ^ . , ? c n (5.22) 
1 p i 'L- 1 ZL- 1 Pi = Ti i p i ,for / = 0, 1. 
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5. Second class of exact solutions 



Since the operator Z is positive, we know it exists an operator Y such that Z = YYK We can 
insert it in p,E _1 ZZ _1 p ; = T] ; p ; and find that it exists W ; , a unitary transformation such that 



W^Z^Pi = y/ffiVPh * = <U- ( 5 - 23 ) 

Moreover, E is full rank. As a result we can decompose Z as Z = PqHT^ZTT 1 po + 
poE -1 ZE -1 pi + piL -1 ZS -1 po + piS-^-Vi- This directly yields 



J]oPo + mPi + V^oli VPWoWi VpT+ VHoi?i VpW/WbVPo (5.24) 



We finally read off Y * as 

yt = ^W^^+ VTiVpi (5-25) 

where W f = w/Wo. 

We now make use of the relation ZE-> = which is equivalent to Y^E? = since AA^ = 

A = for any matrix A. We can explicitly write Y^E? = with = ^/f^oW^ y/po + y/Viy/Pi 
and W = WqW\ . This leads to the statement 

WW 1 " = 1, 

3Y,W such that { YY^ = Z, (5.26) 



In fact, this relation — y^V^y/poE? = v /r/7 v /pT£'? is only possible when — W is a unitary trans- 
formation coming from a polar decomposition of ^/Poy/pT otherwise theorem 13 in chapter 4 is 
violated. Indeed theorem 13 tells us that the product between <2o and Q\ is lower bounded as 

QoGi > r?or?iF 2 (5.27) 
where the equality holds if and only if a unitary operator V arising from a polar decomposition 



VPoVPi = V VPoPi VPo V (5.28) 

satisfies 



V t v ^ v /£7= av^Tv^ (5.29) 

for some a E R + . Moreover F = maxj/ \Tr(U* Jp~Q-Jp~\) | is reached only for unitaries U coming 
from a polar decomposition of ^/po^/p~l■ For any unitary V which does not come from a polar 
decomposition, we then have the strict inequality F > \Tr(V^ ^/po^/pi) |. In other words, if V 
does not come from a polar decomposition then 

W?i^ 2 >^om|Tr(V f VpoVpT)l- (5-30) 



5.2. Maximum rank and a priori probabilities 
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Moreover, since y/po^/E? = cty/piy/E?, the Cauchy-Schwarz (in)equality tells us that 

GoQi = r]QmTr{E,p Q )Tr{E lPl ) (5.31) 

= TfoTh|7>(vVPoE?>/Pi)l ( 5 - 32 ) 
= ^oml^K^VPoVPDI- (5-33) 



Consequently r\QT)\F 2 > <2o£h and the theorem 13 is violated. This implies that —W comes 
from a polar decomposition of y/poy/Pi- At that point, we simply use the equivalence derived in 
chapter 4 

- y/VoW J VPoE? = VmVpiE 7 & { V^l (5.34) 



Indeed Theorem 13 tells us that, for any — W coming from a polar decomposition of \fp~0\fpl, 

- VVoW^y/PoE? = VTiVpiE? & Q°P l = 2^rmF. (5.35) 
And Theorem 16 says that, for any — W coming from a polar decomposition of y/poy/pi, 



p°-^M Fo -° 



Q°P l = l^WViF & { V nF ( 5 - 36 ) 



Pi " JfFi > 0- 



This completes the proof. 



There are at least three consequences to the theorem above. First, it indicates that an 
optimal POVM is, in general, unlikely to have its elements Eq and E\ of maximum rank. This 
comes from the fact that the positivity of two operators po — y^^b an d Pi _ \f^\^ 1 * s on ^ 

possible the middle regime defined by Tr (^ p °) < < Tr( -^ )Pl ) • Second, we can use Theorem 
18 to investigate further the spectrum of an optimal USDM. Last but not least, we can derive 
a new class of exact solutions for the problem of unambiguously discriminating two mixed states. 



5.2 Maximum rank and a priori probabilities 

Theorem 1 8 can be rephrased as 



( p - F^Fq > f rank ( E o) < dim(y pQ ] 

If I V *to is v j i atec i tnen ) or 

Pl ~ V m Fl ~ 1 rank{Ei) < dim(y pi ) 



(5.37) 
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5. Second class of exact solutions 



In this section, we discuss why Theorem 18 suggests that Eq and E\ have maximum rank only 
in a small regime of the ratio between the two a priori probabilities around 1 . 

We already know that the positivity conditions in (5.37) are quite restrictive since they are 
reachable only in the middle regime of the ratio y^- Indeed we repeat here that 

(5.38) 

where Q opt = 2^/r)Q~f\\F is an overall lower bound on the failure probability that cannot be 
reached in the two outer regimes. 

Second, the boundaries of this middle regime can actually be made tighter. Indeed the three 
regimes of the ratio where built considering some constraints on Qq and Q\. Stronger 
constraints means tighter boundaries and the constraints on Qq and Q\ could in principle be 
made stronger if more knowledge on the two density matrices po and pi is provided. 

Let us give such an example of stronger constraints on Qq for, say, a POVM having the 
symmetry E\ = UEqU where U is a unitary transformatiorQ. 

Since E C Jt pi , there exists R > in X Px such that E\ + E? — Pi + R. Moreover the POVM 
element £? is invariant under U since UE?U = U(t — Eq — E\)U = (1 — E\ —Eq) = E?. Hence, 
Eq + E? = U (Ei +E?)U = Po + URU. We therefore obtain the trace equality 

Tr(£ ? ) =2Tt(R). (5.39) 

Indeed Tr(E 1 + E?) = Tr(Pi) + Tr(R) and Tr(£ + E?) = Tr(P ) + Tr(/?) so that Tr(l) + 
Tr(E?) = Tr(Po) + Tr(Pi) + 2Tr(i?). And, for a standard USD problem, the equality Tr(l) = 
Tv(Pq) +Tr(P l ) holds. 

We can now consider Qq. Ei + E? = P\ + R and Tr(.E'iPo) = 0, we can consequently write 

Qq = 77oTr(£ ? p ) (5.40) 
= r]oTr( J B ? po) + T]oTr( J E 1 po) (5.41) 
= r ? oTr(P 1 p ) + T]oTr( J Rpo). (5.42) 

The operator P^PqP^- is a positive semi-definite operator so that its eigenvalues are all positive 
or equal to 0. We can here introduce its smallest non vanishing eigenvalue. It follows that 



We will see in the next section that such a symmetry is possible for USD of two geometrically uniform states. 



5.3. A fourth, incomplete, reduction theorem 
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Go > T?oTr(Pip ) + ri Tr(R)?i min . Together with Eqn. (l539l) this yields 



Q > 7 7o Tr(P 1 po) + ^^Tr( J E ? ) (5.43) 
> T]oTr( J P 1 po) + ^^Tr( J B ? p ). (5.44) 



In other words, for any USD POVM such that E\ = UEqU where U is a unitary transformation, 



Q«>f^ (5-45) 



where A m ,„ = min{Spec(P^- poPj 1 )} . It becomes clear that with more knowledge on the mixed 
states po and pi, we could make the boundaries of the middle regime tighter. The extreme case 
would be a middle regime reduced to y^j- = 1 • These considerations might indicate that, in 
general, Eq and E\ have maximum rank only for some range of the ratio between the a priori 
probabilities around T]i = t\q = 1/2. 



5.3 A fourth, incomplete, reduction theorem 

In the case where y po PI rY Pl = {0}, the maximum rank of the USD POVM elements Ej, i = 0, 1 
is r,-, the rank of the mixed states p ; -. Moreover if not only y po fl 5^ Px = {0} but also Jt Po PI y pi = 
{0} and Jtf Pl fl y pa = {0} then po and p\ have the same rank r in a 2r-dimensional Hilbert space 
and we end up with 

r™ x = r,i = 0,1,? (5.46) 

One can actually use Theorem 18 to study the spectrum of the elements of an optimal USDM. 
In fact, we can state that, for a standard USD problem, if po — y^^b and pi — \J~^F\ are not 
positive semi-definite then the optimal measurement is such that E? possesses one eigenvalue 
equal to 1 and Eq or E\ too. Let us make this result precise in the following theorem. 
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5. Second class of exact solutions 



Corollary 4 A fourth, incomplete, reduction Theorem 

Consider a standard USD problem defined by two density matrices po and pi and their 
respective a priori probabilities rjo and r\\ (any USD problem of two density matrices 
can be reduced to such a form according to Chapter 3). Consider also an optimal 
measurement {E ? 1 ,E° pt ,E° pt } to that problem. Let Fq and F\ be the two operators 
\/^/poPu/po and \J ^/pipoy/pi. The fidelity F of the two states po and pi is then given by 
F = Tr{F ) = Tr(F l ). 




\e) G y po and \ e') G J^ such that 



rjopt 
C 



e) 



E opl 



ET\e) = \e) 
e>) = \e>) 
e)~- 



E opt | 



E opt \e') = 0, 



or 




First let us note that this theorem makes this assumption of a standard USD problem. It is in 
principle not necessary to make such an assumption to derive the existence of some eigenvector of 
E?, Eq or E\ with eigenvalue 1 since Theorem 18 is valid for any pair of density matrices without 
overlapping supports. Nevertheless, this theorem aims to be a 'fourth' reduction theorem. It 
means in particular that, for any given USD problem of two density matrices, we would like to 
apply our 'four' reduction theorems and always end up with the optimal USD measurement. 

The above theorem is a kind of incomplete reduction theorem. A reduction theorem is a 
theorem that allows us to decrease the size of the USD problem by splitting off some subspace 
onto which no optimization is needed. To have a complete reduction theorem here, we would 
need to characterize \e) and \e') without solving the whole optimization problem. But only 
the existence of \e) and \e') is so far ensured. If such a reduction theorem were found then 
we would have a recipe to solve any USD problem. Let us assume that \e) and \e') are fully 
characterized and let us start from a general USD of two mixed states. We use the three first 
reduction theorems to make it standard. We then check whether the two operators po — y^^b 

and pi — J j^-Fi are positive semi-definite. If yes then we know the optimal failure probability 



5.3. A fourth, incomplete, reduction theorem 
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as well as the optimal measurement to perform since this case falls into the first class of exact 
solutions (middle regime). If the two operators po — y^^b an d Pi — yrj^i are not positive 
semi-definite, we can use our last reduction theorem to get rid of two dimensions. At that 

point, we check again the positivity of the two operators p' Q — y^f^o an ^ P\ ~~ v/^f^i °^ 
the reduced problem. We see here a constructive way to solve any USD problem. If the two 
operators p' Q — \ j^~Fo an d p[ — \ never happen to be positive, we end up with only two 

pure states and can finally find the optimal measurement (see Fig. l5.ll) . The only problem in that 
nice picture is that we only know that \e) and \e') exist but we cannot until now characterize them. 



Standard form 



rr j 



(r) (r) 
(r) (r) 

pi , m 



Fi r) > 



— 1> > n 

i) Ml — u 



Pi 



> 



Q ' = 2ym m F 



(2) 
Pb 



P? 



R*l F m 
V 



Q (2) - z^W^ 



■laeger and SMmony's 
solution for 
two pure states 



Figure 5.1: A constructive way to solve any USD problem (the exponent ^ denotes the rank of 
the density matrices after reduction) 



Here comes another important remark. There are only two ways to find a complete charac- 
terization of the two eigenvectors \e) and \e'). The first is to consider a low dimensional USD 
problem. The second is to consider a highly symmetric problem. The former case simply is 
the two pure states case. Indeed, either the operators po — y^^O an d Pi — y^^i 316 positive 
semi-definite or we have \e) G ^p 0/1 and \e') G ^p 0/I , eigenvectors of E? and £i/o- In only two 
dimensions, there is no freedom and \e) and \e') must be |Yq/i) and | Y^y r ) . If we are interested 
in higher dimensions, we use some symmetry to give us enough constraint to fully characterize 
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\e) and \e'), we can go up to four dimensions. This is the object of our last section. Before that 
let us prove Corollary 4. 



Proof of Corollary 4 To prove this corollary, we begin with the statement given in Theorem 
18 for two density matrices po and p\ with same rank n in a 2n-dimensional Hilbert space. The 
maximum rank of Eq and E\ then equal n. Let us for example consider that rank(Eo) < n. The 
other option corresponding to rank(E\) < n follows the same argumentation. Because of the 
completeness relation E? + E\ + Eq = 1 fulfilled by the POVM elements, we have, onto the 
subspace yp , the following equality PqEiPq -\-PqE\Pq -\- PqEqPq = Pq. However, y% x G J^pt so 
that we are left with 



P E ? Po + PoE Po = P . (5.48) 

Furthermore, in PqEqPq's eigenbasis, we have PqEqPq = YJiZ] M\^i){^d\ since Eq is of rank 
n — 1 and Pq = E^=i 1^') (^'1 + \ e ) i e \ where \e) completes the n dimensional orthogonal basis of 
yp . As a result, E?\e) = (1 — Eq — Ei)\e) = \e) — — and \e) is an eigenvector of E? with 
eigenvalue 1. 



We can actually go one step further. Since the completeness relation is already fulfilled 
onto the subspace spanned by \e) and \e'), no optimization is required onto it and we can split 
it off from the original USD problem. The remaining USD problem to optimize concerns 
Pq and p[ originated respectively from the density matrix po and pi. Moreover, p' has rank 
n — 1 while p[ has rank n. We can indeed denote by the subspace spanned by \e). The 
reduced Hilbert space is 1 5?\ e } and y po , the support of po, looses one dimension. Thanks to 
the second reduction theorem, we can reduce this problem to the one of two density matrices 
of rank n — 1 in a Hilbert space of dimension In — 2. Indeed, the subspace J^ p i n is 
one dimensional and leads to the detection of p[ with unit probability. We call \e') the unit 
vector spanning this 1 -dimensional subspace. We are left with a reduce USD problem in a 
In — 2 dimensional Hilbert space. Importantly, \e') is in n =5^ C =5^7 = Indeed, 

w = y Po © ^pi = y P ' Q © y\ e ) © y p \ ™ that, in m" = je/y\ e) , ^ = y^. 

In other words, if po — y^^O an d Pi — \J^[F\ are not positive then it exists \e) in S^p^, 
eigenvector of E? with eigenvalue 1 and \e') in ^p , eigenvector of E\ with eigenvalue 1. 
Without the assumption that rk(Eo) < rk(pj), we have in general that if po — y^^b an d 

Pi — \J^F\ are not positive then there exists \e) in either yp or yp v eigenvector of E? with 
eigenvalue 1 and \e') in either J(f P{) eigenvector of E\ with eigenvalue 1 or JXf Pl , eigenvector of 
Eq with eigenvalue 1 . The completes the proof. ■ 
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The third consequence of Theorem 18 is the derivation of the optimal USD measurement for 
any pair of two geometrically uniform states in four dimensions. 



5.4 Second class of exact solutions 

Geometrically uniform states, or GU states, are a generalization of symmetric states ll50l I5TI l52l 
|22][53l[36]|. While symmetric state are generated from one generator state and a single unitary 
transformation, GU states are generated from one generator and a group of unitaries. They are 
interesting for both practical and theoretical considerations. On the practical side, real applica- 
tions often exhibit strong symmetries like GU symmetry^. On the theoretical side, this symmetry 
allows us to seek for simpler conditions and then new results. Actually Eldar proved that the op- 
timal measurement to unambiguously discriminate geometrically uniform states can be chosen 
geometrically uniform, too. This result allows us to derive now the general solution for unam- 
biguously discriminating any pair of GU states in four dimension. Next we give the mathematical 
definition of the geometrically uniform states before presenting the optimal failure probability for 
unambiguously discriminating two geometrically uniform states in four dimensions and the cor- 
responding optimal measurement. 



5.4.1 Geometrically uniform states 

A set of GU state is a set of mixed states {pi}, i = l,...,n such that p,- = £/,p£/ ( where p is 
an arbitrary density matrix called the generator and the set {£/,}, i = \,...,n is a set of unitary 
matrices that form an abelian group. In order not to break the symmetry of the states, we assume 
that all their a priori probabilities are equal to -. 

A consequence of the group structure of the set {£/,} is that we can always consider U\ as the 
identity, and pi as the generator for a given set of GU states. We can therefore always write two 
GU states as po and p\ = U poll where U is an involution (i.e. a unitary transformation U such 
that U 2 = 1) with t\q = r\\ = ^. Let us note that two GU states are two symmetric states since 
only a single unitary is needed. 

In the next section, we give a second class of exact solutions for USD of two generic density 
matrices. We provide the optimal failure probability as well as the optimal USD measurement 
for any two GU states in four dimensions. 



2 In a cryptographic context, the bit value states and basis states in the BB84-type protocol using weak coherent 
pulses and a phase reference exhibit such a GU symmetry. 
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5.4.2 Optimal unambiguous discrimination of two geometrically 
uniform states in four dimensions 

Theorem 19 Optimal unambiguous discrimination of two geometrically uniform states in four 
dimension 

Consider a USD problem defined by two geometrically uniform states po and pi of rank two 
with equal a priori probabilities and spanning a four- dimensional Hilbert space. Let Fq and 
F\ be the two operators a/ ^/Popi^/Po and 1 ^7p7p~o~\7pT- The fidelity F of the two states po 
and pi is then given by F = Tr(Fo) = Tr{F\). We denote by Pq and P\, the projectors onto 
the support of po and p\. The optimal failure probability Q° pt for USD then satisfies 

1. 2° pt = F ifp -F >0 (5.49) 



2 - Q ° Pt ~ 1 {xlPolx) if {7pec(P^UP^) = {a,-b}, a,beR+ 

3. Q ? 1 = 1 otherwise. 

with P^UP l ± = a\0)(0\~b\l)(l\and\x) = -^(g-^^NO)) y^| ) + ^| 1». 

The POVM elements that realize these optimal failure probabilities are given in the different 

cases by 



l.E = Z- l ^(p -Fo)^T,- 1 (5.50) 
Ei = UEqU 
E? = 1-Eq-UEqU 




3. E = 
Ei = 
E? = 1. 



Proof We consider a USD problem defined by two geometrically uniform states po and 
pi = U PqU, U 2 = 1, of rank two, spanning a four-dimensional Hilbert space. This means in 
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particular that J^ Po n y pi = {0} and r'™ x = r™ x = r'™ x = 2. 

Due to the symmetry of the states, we also notice that po — Fq = pi — F\. Note that the a 
priori probabilities are equal in order not to break the symmetry. Moreover, thanks to Eldar ll36ll . 
we can choose the optimal USD measurement to be GU, too. Thus the POVM elements are such 
that 

Eq , (5.51) 
El =UE Q U, 
E? =UE?U. 

The statement in Theorem 16 for equal a priori probability 

2° pt = F # (5-52) 

then reduces to 

Q opt = F & p -F >0. (5.53) 

Note that we are not interested in the equivalence. The implication from the right to the left 
is the only important direction for our purpose here. In that case we need the assumption 
^po n ^pi = {°} t0 prove that: If p - Fq > then Q opt = F. Without this assumption, only the 
other direction is true. 

If po — Fq ^ 0, Theorem 18 tells us that the ranks of the POVM elements Eq and E\ are 
not maximum (Eq and E\ have the same rank because of the symmetry). As a consequence, if 
Po - F ^ then rank(E ) = rank{E\) < 2. It follows that if p - F ^ then the two POVM 
elements Eq and E\ have either rank 1 or rank 0. If rank(Eo) = rank{E\) = then E? = 1 and 
2=1. Let us now focus on the remaining case rank(Eo) = rank{E\) = 1. 

Let us now prove that a measurement with rank(Eo) = rank{E\) = 1 and rank(E?) < 2 is 
necessary a projective measurement with rank(E?) = 2. We can introduce the unit vectors and 
real numbers |jc) G J^ Pl , \y) G Jtp , x and y such that 

E = x\x)(x\,E l =y\y)(y\. (5.54) 

We call S^xy the two dimensional subspace spanned by |jc) and \y), Pxy the projection onto it and 
P^y the projector onto its orthogonal complement. By definition of the subspace S^xy, 



(5.55) 
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Therefore rank(P^,E?P^) = rank(P^) = 2 and Ej must be at least of rank 2. However 
rank(E?) < 2. Therefore rank(E?) = 2 and 

E? = P^. (5.56) 

We can now consider the subspace y xy only. On that subspace, we have 

£ + £i=/> l7 (5.57) 

that is to say Pxy = x\x)(x\ + y\y)(y\. Since P xy is a projector, P^ = P 2 and it follows that 
jc|jc)(jc| +y|v)(y| + xy(y|x)|y)(x| +xy(y|x)|y)(x| = x|x)(x| +y|y)(y|. The off-diagonal terms are 
equal if and only if (y\x) = while the diagonal terms are equal if and only if x = y = 1. The 
POVM then is a projective measurement with rank(E?) = 2. 



We now give the optimal USD measurement for a GU projective measurement. Since the 
measurement is made of projectors, we have Tt(EqEi) = which is nothing but (x\U\x) = 0. 
Because |x) lies in J^ Px , this relation is equivalent to 

(x\Pj-UPj-\x) = 0. (5.58) 

P^UP^- is a Hermitian operator and therefore owns real eigenvalues. Note that if P^UP^- must 
be of rank 2 since U is full rank. Thus we denote a and c the two eigenvalues of P^UP^~ and |0) 
and 1 1) its two eigenvectors. In this eigenbasis, \x) G J^ Pl can be expressed as 

W = (p) (5 - 59) 

This leads to {x\P^UP\\x) = \a\ 2 a + |/3| 2 c. Importantly this scalar product can only vanish if 
a > and c < 0. We call -c = b > such that, in {|0), |1)}, 

PtUPt =(l° b )- (5-60) 

If we include the normalization of \x), we end up with a system of two equations. This system 
simply is 

|a| 2 a+|/3| 2 c = 

l«| 2 +|/3| 2 =l (5 - 61) 
and admits a family of solutions parametrized by a phase <3>: 

{a= t e =,/3 = -^^=}. (5.62) 
1 Vl+a/b Vl+b/a S 
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In the basis {|0), 1 1)} we can therefore write 

\x) = ( ^ l \ a/b ) . (5.63) 

We can use again the fact that we are interested in the optimal measurement. Note that we 
already considered optimality to state than if po — Fq ^ then the POVM is either {Eq — E\ — 
0,E? = 1} or a projective measurement. Indeed Theorem 18 is only concerned with optimal 
USD POVM. So far, \x) is valid for any USD measurement such that Eq = \x){x\, E\ = UEqU 
and E? = t — Eo — UEqU. Let us now find the optimal one. To do so, we evaluate the success 
probability P^ccess- Because of the symmetry of the two GU states, Tr(£bPo) = Tr(£ipi) and 
the success probability Process = ^Tr^opo) + \Tr{E\p{) for unambiguously discriminating the 
two GU state po and pi takes the form 

PTccess = Tr(£ po) = (x\p \x). (5.64) 

After calculation, we obtain 

PTccess = ^(MO|Po|0)+a(l|po|l)+2v / ^i?e((0|p |l)^)). (5.65) 

We choose the phase O to maximize this success probability Process- That is why we choose <I> 
such that /te((0|po|l)e f<I> ) = |(0|p |l)|. Therefore, O must be -Arg((0\po\ 1)) and 

/ e -^g((o|p |i» \ 
|x) = I vTWfc | . (5.66) 



Vl+b/a 



This completes the proof. 



This theorem leads to a fundamental question: 'Is it possible to find a unified expression for 
the failure probability QT In the first class of exact solutions, we can write the three failure 
probabilities of the three regimes as 

Q = ar] F+-ri 1 F 

with the above-mentioned a. But we do not really expect the bounds in the outer regimes to be 
often optimal (see discussion in section 5.2) so that this expression does not seem so fundamental. 
More significatively, for the second class of exact solutions, no unified expression of the failure 
probability exists. In higher dimension (dim(J^) > 4), the number of cases for the optimal 
failure probability Q might become very large. If this is the case, a unified expression for Q 
would be a pre-condition to find the general solution to USD of two density matrices. 

In the next chapter we analyze an application of both theoretical and practical interest. In 
fact, we consider the Bennett and Brassard 1984 protocol (BB84 protocol) implemented through 
weak coherent pulses with strong phase reference. This represents the first solved example of a 
non reducible USD problem. 
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Chapter 6 



Application of the second class of exact 
solutions to the BB84 protocol 



In 1984, Bennett and Brassard proposed a protocol to distribute a unconditional secure private 
key between two parties over a public channel in order to allow a secure communication. This 
proposed Quantum Key Distribution protocol, the so-called Bennett-Brassard 1984 (or shortly 
BB84) is here unconditional secure because of the laws of nature (quantum mechanics) and 
not anymore because of the assumption of a limited computational power of some hypothetical 
eavesdropper. In the standard BB84 protocol, Alice sends one of the four states {0, 1,+,—} to 
Bob. Here {0, 1} and {+,— } are orthogonal pairs and and + correspond to the bit value 
while 1 and — correspond to the bit value 1 . Bob then detects the signal sent in one of the two 
bases {0, 1} or {+,—}. 

In this thesis, we consider the implementation of a BB84-type protocol that uses weak co- 
herent pulses with a phase reference. In that scenario, Alice sends one of the four states 

The bit value is encoded in the sign of the coherent states that is to 

say |-^=) and |^=) correspond to the bit value 0, |^|) and \ correspond to the bit value 1. 
Moreover the phase i plays the role of the basis in the standard BB84 protocol. Firstly let us note 
that the factor 4= in the amplitude comes from the technique used to implement the polarized 
coherent states. Secondly the first mode is common to the four signal states. This mode 
is therefore irrelevant for the following analyze. Furthermore it is worth noticing that the states 
corresponding to the bit value and 1 are not orthogonal since 



(6.1) 
(6.2) 
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This QKD protocol is therefore not the standard BB84 protocol. It remains that two important 
question can be addressed. 



With what probability can an eavesdropper unambiguously distinguish the basis of the 
signal? 



With what probability can an eavesdropper unambiguously determine which bit value is sent 
without being interested in the knowledge of the basis? 



In fact the first question refers to the unambiguous discrimination of the two basis )} 
and {|yf>}- Therefore we can build a mixed state po that corresponds to the basis {|^r)} an d 
a mixed state p\ for the basis {1^^)}. We end up with 



Po 
Pi 



1 

2 

1 / 

2 



a a -a -a 

7! v! 7f )( 7f 



V2 



(6.3) 
(6.4) 



where we ignore the irrelevant first mode. 



The second question refers to the unambiguous discrimination of the two bit value mixed 
states. We can for that case build the two density matrices 



Po 
Pi 



1 

2 

1 / 

2 



a a . ,ia,.ia 

71 VI 75 71 



a. .— i 



y/2 W2 y/2 y/2 



(6.5) 
(6.6) 



where we again ignore the irrelevant first mode. 



The states )}> {1^7^)1 are f° ur linearly independent pure states. Therefore they span a 
four dimension Hilbert space. In the next section we will express the four density matrices above 
in that four dimensional Hilbert space and prove that they are GU states. After that, we will solve 
the two USD problems arising from the two questions mentioned. It turns out that the first case 
is reducible to some pure state case while the second one requires our last theorem to be solved. 
Let us now start with the explicit expression of these four mixed states. 
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6.1 Two geometrically uniform states in a four- 
dimensional Hilbert space 

A coherent state of amplitude a can be written as a poisson distribution of photon number in the 
polarization mode cfi as 

where |0) denotes the vacuum state. Moreover, the four signal states | ± a), |/± a) are coherent 
states in four different polarizations: ±45° and circular left or right. These polarizations are 
expressed in terms of two orthogonal polarizations b[ and b\ as 

4 = -^p\+ b i)> ( 6 - g ) 

a\ = ±={b\ + ib\), (6.9) 
4 = -^(b\-bl), (6.10) 



V2 
1 



a\ = J-(b\-ib\). (6.11) 



|Yo) = 


I—) 




|Yi> = 






I—) 


TO = 


I—) 




|Y ) = 


I—) 





Consequently, we can write the four states as 

(6.12) 
(6.13) 
(6.14) 
(6.15) 

The first mode is common to the four states and therefore will be left out. In the phase 
space, these four states are generated from |Tn) and a rotation of angle j. This means they are 
symmetric states and we can write them in a suitable basis following Chefles et al. [T9|. The 
idea is that n symmetric states can always be written in an orthonormal basis { } as 

|Y*>= I*;>- (6-16) 

7=0 

Note that the phase of the complex numbers Cj is not relevant since we can absorb it in the defi- 
nition of the basis elements Actually the modulus of the coefficients cjs can be expressed 
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lfT9l as 



This leads in our case to 



-2i% 



j(k-k') 



Lid 



\co\ 


1 

— ~l= e 


M 


lcosh(^) 


+ COS' 


Us 


ki 


1 

~V2 e 


_M 


Jsinh(^ 


) + sin< 


,/x 


\ c 2\ 


1 

= vT 


M 




— COS' 


,jU 

-2 j 



c 3 



-J=e 4^jsinh(^)-sin(^). 



(6.17) 

(6.18) 
(6.19) 
(6.20) 
(6.21) 



(6.22) 



where /i = |a| 2 stands for the mean photon number. Moreover, in the basis {|<E>j)}, the unitary 
transformation acting on | Yo) that generates the other three states is 

/ 1 \ 
0/00 
0-10 
\ — i J 

such that K 4 = 1. The four symmetric states (see Fig. 16. II) are then expressed as 

/c \ 
c\ 

\C 3 ) 



l*o> 



(6.23) 



|Yi> = *|Yo> 



|Y 2 ) = £ 2 |Y > 



|Y 3 ) = £ 3 |Y > 



/ co \ 

ic\ 

-C2 
V " /C 3 / 
/ Co \ 

-ci 

C2 

V " C 3 / 
/ CO \ 

—ic\ 

-C2 

V /C 3 / 



(6.24) 



(6.25) 



(6.26) 
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At that point, we are ready to write the four density matrices corresponding to the basis mixed 
states and the bit value mixed states. 
The basis mixed states (see Fig. I6.2I) 



po = i(|Yo)<To| + |Y 2 )<Y2|), 
p, = i(|Y 1 )(Y 1 | + |Y 3 )(Y 3 |) 



(6.27) 
(6.28) 



are by construction of rank 2. 

They can be written in a four dimensional Hilbert space spanned by the four linearly inde- 
pendent states |Y,), i = 0, 1,2,3 as 



Po 



and 



Pi 





/ c 2 





C0C2 


\ 











c l 













cqc 2 
















\ 


ClC 3 











/ 







-c c 2 







\ 







r 2 
c 1 





-C1C3 






cqc 2 





r 2 









V 





-C1C3 





r 2 




/ 



(6.29) 



(6.30) 



where we choose all the coefficients c,- to be real. 
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Ol*3> 



Figure 6.2: Pairing of the four symmetric states for the basis mixed states 



Thanks to Eqn. (6.27) and Eqn. (6.28), we clearly see that 

Pl = Kp tf = tfpoK. 
Moreover, we can calculate that KpoK = K^poK^ in the following calculation. 

Kp K = ^(K\Y )(Y \K + K\Y 2 )(Y2\K) 

= ^^\Y 2 )(Y 2 \KUKl\Y )(>¥o\K< 
= K^p K\ 

The consequence is that we can construct two new unitary matrices which are involution^ such 
that pi = U±PqU±. This two involutions are given by 



(6.31) 

(6.32) 

(6.33) 

(6.34) 
(6.35) 



U± = — z — ±i — - — 



(6.36) 
(6.37) 



We can choose to use in the following calculation 

/ 1 \ 
10 
0-10 
\ -I J 

l A unitary transformation U is called an involution if and only if U 2 = 1. 



u = u 



(6.38) 
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We have finally written the basis mixed states as po and p\ = U PqU where U 2 = 1. This means 
that the question 'With what probability can an eavesdropper unambiguously distinguish the 
basis of the signal?' is related to the unambiguous discrimination of two geometrically uniform 
mixed states in dimension four. The choice of such a involution matrix will simplify the next 
calculations. Finally, in the four dimensional Hilbert space, we see that 



such that rank(po + Pi) = 4 = rank(po) + rank(pi). The two GU states po and p\ do not have 
overlapping supports and we can apply Theorem 19 about USD of such a pair of states. The bit 
value mixed states (see Fig. 16.31) are also rank two matrices by construction. They can be written 

as 



Po + Pi 



( c\ \ 

c\ 

c\ 

V c\ ) 



(6.39) 



Po 



Pi 



i(|Y )<Y | + |Yi)<Yi|), 
i(|T2><T 2 | + |T3>(T 3 |). 



(6.40) 



(6.41) 



Im(a) 




^Re(a) 




Figure 6.3: Pairing of the four symmetric states for the bit value mixed states 
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In terms of the coefficients c ; 's, we obtain the following form i 
space spanned by the states |Y,), i = 0, 1,2,3: 

/ c 2 Q ^cqCi if c c 



in the four dimensional Hilbert 



Po = 



if CiC 





1-i 



V V C 3 C 



c 2 
^C 3 C 2 





if c 2 c 3 

J- 



(6.42) 



/ 



and 



Pi = 



/ r 2 


-if cic 




■V c o c i 

C 2 ~ "f-ClO 

— C2C1 c 2 



C C 3 \ 



-^f C 2 C 3 

/ 



V -V C 3 C — f C 3 C 2 C 2 , 

It is unfortunately impossible to choose the phase of the coefficient c,- 
matrices. Therefore we simply choose all the coefficient c, to be real a 



(6.43) 



2 -3 / 

hase of the coefficient c,- so that po and p\ are real 
coefficient c, to be real and we end up with 



Po 



if ClCQ 



y if c 3 c 









I 2 i ClC 2 

r 2 

ifc 3 C 2 



^ C C 3 \ 



if C 2 C 3 

„2 



(6.44) 



and 



Pi 



/ 



V 



p 

if C\CQ 



^ C3C0 



-if C Cl 

r 2 
c l 

"if C 2 C1 







l-( 



-if c c 3 



cic 2 



i fc 3 C 2 





I-; 



QC3 



(6.45) 



The involution connected po and pi simply is 

/ 1 



K 1 



\ 



0-100 



1 









1 



(6.46) 



/ 




\0 

Of course, the question 'With what probability can an eavesdropper unambiguously determine 
which bit value is sent without being interested in the knowledge of the basis?' is also related to 
the unambiguous discrimination of two geometrically uniform mixed states in dimension four. 
Here again, the sum po + Pi in the four dimensional Hilbert space is given by 



f cl \ 



P0 + P1 



d 



.a 
1 





(6.47) 



V 4 ) 



6.2. USD of the basis mixed states 
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implying that the two GU states po and p\ do not have overlapping supports. Consequently 
Theorem 19 can be used. 

Actually one could consider a third USD problem coming from the pairing of the four states 
Yj (see Fig. 16.41) . This last case is concerned with the unambiguous discrimination of the two 
mixed states p = \ (|Y )(Yo| + |Y 3 )(Y 3 |) and pi = \ (|Yi)(Yi| + |Y 2 )(Y 2 |) but this case is 
similai@ to the previous case. Indeed one can go from the former to the later case by using the 
unitary K 2 . This is not the case between the two problems of unambiguously discriminating the 
basis states and the bit value states. 




Ol*3> 

Figure 6.4: Third possible pairing of the four symmetric states 



6.2 USD of the basis mixed states 

Let us repeat that the two density matrices to unambiguously discriminate are 



Po 



/ 


r 2 





cqc 2 





\ 







r 2 
c l 





cic 3 






CQC2 





r 2 
6 2 







\ 





ClC 3 





r 2 


/ 



(6.48) 



2 unitary equivalent 
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and 



pi = Up U 



( 


c 2 





-c c 2 





\ 







c l 





-C1C3 






~CqC 2 





c 2 

6 2 







\ 





~C\C3 





r 2 
c 3 


/ 



(6.49) 



with 



U 



( 


1 











\ 







1 


















-1 







{ 











-1 


) 



(6.50) 



With a bit of concentration, one can realize that these two density matrices are block diagonal. 
Indeed, we can use the permutation matrix 



/ 1 








1 





o\ 






(6.51) 



and obtain 



/ 


c 2 


±coc 2 








\ 




±C C 2 


c 2 


















c 2 


±ClC 3 




V 








±ClC 3 


r 2 
c 3 


/ 



(6.52) 



This already tells us that we can analytically solve this problem which is reducible to some pure 
states case. Indeed po and pi are block diagonal where each block is two dimensional. We 
will nevertheless use the non reduced density matrices to find the optimal USD measurement. 
The reason is that, as we will in the next paragraph, we can compute the operator po — Fq and 
check its positivity for any value of the amplitude a. Note here that the spectra of po — Fq and 
Pi — Fi are identical since p\ — F\ = po — Fq for two GU states. With that, we have the optimal 
failure probability as soon as the optimal measurement. Again, we could use the second and 
third reduction theorems but the present example gives us the opportunity to use other tools. 
We now focus our attention onto po only since p\ is similar to it. The density matrix 



/ 



PpoP 



C C 2 






coc 2 
c 2 









c 2 
c l 

ClC 3 






c\c 3 



\ 



(6.53) 
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can be easily diagonalized using the block diagonal unitary matrices 



/ 



pu q p 



CO 



V 



V /c ol 2+c o 






C2 



V c o+ c o 









c\ 



g3 






C3 



(6.54) 



If is not too difficult to find that the eigenvalues of PpoP are therefore given by 



4+4 
4+4 



which gives, in terms of the mean photon number /i 



M),l 



(6.55) 
(6.56) 



(6.57) 



If we undo everywhere the permutation matrix P, the density matrices po,i can obviously be 
diagonalized with the help of the unitary transformation 

/ 



U 



V c o+ c o 








C2 



C2 



V 



V c 0+ C 





and its square root takes the form 

/ 



V c 0+ C 



CQ£2 



\ 



o 










_£l£3 



V c 0+ c 



-go 





C3 











(6.58) 



CQC 2 



V c 0+ C 





v 7 ^ 



V C 0+ C 





(6.59) 



The next step is to calculate the operator Fq 



through the relation ^fpi = U y^poU. As a result, the equality ^/pQ^/pi 

y/frU^ = F VU. 
In the po's eigenbasis, we obtain 

t/ovWvWo = UqFqVUUq = UqFqUqT 



PoPiVPo- ® m two GU states are related 
FqV leads to 

(6.60) 



(6.61) 
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where T = UqVUUq is a unitary transformation. One can calculate the operator Uq^/PqU ^/PqUq 
and find 



/ c 2 Q -c 










\ 



c\-c\ 




V 



(6.62) 



/ 





which is always positive if multiplied by some signature matrix 

/ ±1 \ 

±1 







(6.63) 



Note here that, in terms of the mean photon number /I, the quantities Cq — c| = e^~cos^ and 
c\ — C3 = e^sirij are not always positive. In the end, the positive operator Fq is of the form 



/ \cn — c 



UqFUq 



\ 



V 









4\ 









j 



(6.64) 



The explicit form of the unitary V is only relevant to calculate the elements of the optUSDM. 
But our first goal is to find the spectrum of the operator po — Fq. For that, four cases are to take 
into account depending on the sign of Cq — c\ and c\ — c\. 



Everything is gathered to obtain the explicit form the operator po — Fq in the eigenbasis of 
Po- Indeed, we have 



U (po-F )Uo 



( cl + cj \ 
c\ + cl 



/ \c 2 -c 2 \ 



+ 



J 



V 



000 
000 

^ max{c^,c\} \ 

max{c\,c%} 














>0. 



\ 





J 



(6.65) 



(6.66) 



The spectrum of the operator po — F0 is positive for any value of the mean photon number ji. 
As a consequence, the optimal failure probability Q reaches the lower bounds F = Tr(Fo) = 
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|cq — c\\ + \c\ — eg]. In terms of the mean photon number /1 (see Fig. I6.5I) . the optimal failure 
probability is given by 



2 = e 2 ^|cc»5— | + | sin — 



(6.67) 




Figure 6.5: Optimal failure probability for USD of the basis mixed states 
Let us note here that if we were interested in the unambiguous discrimination of 



Po 
and pi 



1 

2 
1 



\a){a\ + 1 - a){-a\) 
\ia)(ia\ + I - ia){—ia\) . 



then we would find 



Q = e V- (\cosn\ + \sinn\) . 



(6.68) 
(6.69) 

(6.70) 



Let us conclude this section and this example by adding that we can give the optimal mea- 
surement to achieve Q = F. Indeed, the useful matrix E is diagonal and therefore its inverse 
simply is 



'-1 



/ 


Co 2 











































V 













) 



(6.71) 



In the four different cases parametrized by the signature T , the elements of the optimal POVM 
are finally given by 

E = E- 1 v ^(po-F ) v ^E- 1 (6.72) 
Ei = UE U (6.73) 
E ? = E- 1 ( v ^+ v / PTV t )^o(v / P^+^v / PT)^ 1 (6-74) 
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where all the different matrices involved in these equations are perfectly known. This concludes 
this section and the first example. 

6.3 USD of the bit value mixed states 

The second case corresponds to the unambiguous discrimination of the two density matrices 



Po 



/ c 2 




^3 C 



c 2 
^C 2 d 







c 2 







(6.75) 



and 



pi = Up U 



cic 



V C 3 C 



C 2 

-^C 2 C! 







- I y i ciC2 
r 2 



CQC 3 \ 



„2 



(6.76) 



with 



U = K 2 



( 1 







-1 









1 











\ 



1 



(6.77) 



/ 



This USD task is far more complicated than the first one. It is difficult to find the unitary 
transformations to diagonalize po and p\ and therefore the square root of those states as well 
as Fq and F{ cannot be easily expressed. We have to resort to a particular decomposition of the 
two states po and p\. This decomposition allows us to diagonalize the operator po — Fq in an 
unknown basis and find its spectrum. First we review some relevant properties of the density 
matrices po and pi. Next, we solve the unambiguous discrimination of these two GU states. 

Actually one can write 



p =APA 

where A is a real diagonal matrix and P = ^- a pseudo projector. They are defined as 

/ c \ 

ci 

c 2 

V c 3 j 



(6.78) 



(6.79) 
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and 



/ 


1 


1-2 
2 





1+i 
2 


\ 




1+i 
2 


1 


l-i 
2 












1+i 
2 


1 


l-i 
2 




V 


l-i 
2 





1+i 
2 


1 


/ 



(6.80) 



Here come three remarks arising from this decomposition. First of all, let us note that they 
commute since they are both diagonal. Due to the symmetry between po and pi and to the 
commutation between A and U, we have pi =UAPAU = AUPUA. Second of all, we can consider 
the sum of the two density matrices po and pi. We have po + Pi = APA + AUPUA = A(P + 
UPU)A and P + UPU = 21. Thus 

p + Pi=2A 2 (6.81) 

and we could denote A = y^. The last remark is the more important. Actually Tr(P) = 4. This 
is not a lot but it implies that P is equal to twice a two-dimensional projector. As a matter of fact, 
there exists a unitary transformation W so that 



WPW 



t 



/2 












o\ 




2 / 



(6.82) 



Such a unitary matrix can be given by the Discrete Fourier Transform 

/ 1 1 1 1 \ 



W 



\ 



1 



1 1 



1 



(6.83) 



/ 



The interest of the decomposition provide in Eqn.(6.78) is that it allows us to write po = APA 
in an unknown basis better suited to investigate the spectrum of the operator po — Fq. Indeed we 
can write 



Po 



APA 
AP PA 



PoRqRoVPo 



(6.84) 
(6.85) 

(6.86) 



where we introduce the unitary transformation Ro such that ^ = Ro^fpQ. Consequently, we 
obtain 

„ „ t PA AP 
RoPoRl 



V2V2 
PA 2 P 



(6.87) 
(6.88) 
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For the unambiguous discrimination of the two basis mixed states, we knew the unitary trans- 
formation Uq that diagonalizes po- It was possible to write Fq and finally express the operator 
Po — Fq. Here we can not directly work with the eigenbasis of po- Instead, we try to use the ma- 
trix RqPqRq, knowing only the existence of this unitary transformation Rq. We are only interested 
in the spectrum of po — Fq and the precise form of Rq is finally irrelevant as long as it permits 
us to find the spectrum of po — Fq. Nevertheless, we must say, that the explicit expression of the 
POVM elements will not be provided since, it that case, we do need to know Rq. Moreover, as 
we will soon see, we will not be able to calculate the complete expression of Q for all the regime 
of the mean photon number /i. 

Let us now calculate the spectrum of po — Fq. We first apply the Fourier Transform W onto 
RqPqRq to end up with 



WRqPqRIw^ 



l\VPA 2 PW^ 



( cl + cj + cl + cj c\ + ic\-c\-ic\ \ 
















(6.89) 



(6.90) 



\Cq- icj -cl + ic\ eg + c \ + c\ + c\ J 



Actually, since the states po is normalized, we have c\ + c\ + c\ + c\ = 1 and therefore 



WRqPqR^ 



( 1 A \ 





A* 1 



(6.91) 



where 



A 



'A- 



4) 



+ i(cl 



4) 



(6.92) 



In fact, a Hermitian matrix of the form 



a 

be'* 



be'* 
a 



(6.93) 



with a, b and real and positive, has for eigenvalues 

A± = a ± b 

and for eigenvectors 

|v±) 



V2 V 1 



(6.94) 



(6.95) 
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Here we are only interested in the spectrum of pQ. The formula above gives us its eigenvalues as 

1±IM 



a± = (6.96) 

= — . (6.97) 

As for the basis mixed states case where we calculate the operator Uq^/PqU \/PqUq, we now 
consider the operator WRq^/PqU y/poR^W^ . Actually this operator is of a similar form than 
WRqPqRIw^. Indeed we obtain 



/ ( c 2 + c 2 )-(c 2 + c 2 ) 



WRo^U^Rlw 1 



\ 






-A 









-A* 





\ 



[c\ + c 2 ] 



(6.98) 



Sl + c 2 ) ) 



Thanks to Eqn.(6.94), we find that its eigenvalues are 

y ± = (c 2 + c 2 )-(c 2 + c 2 )±\A\. 



(6.99) 



Moreover, with the help of Eqn.(6.95), we obtain the unitary that diagonalizes the operator 
WR 0y /foU ^/p^W^ . This unitary is of form 



/ 



V2 



1 







A* \ 





(6.100) 



A* 

"W 

V2 

y/2 

\ 1 1/ 

If we replace the coefficients c ; by their expressions in term of the mean photon number /i, we 
end up with 

-JL 



( 



\ 

























(6.101) 



The eigenvalues in the top left corner is always positive while the eigenvalue in the bottom right 
corner is always negative. Therefore the operator Fq in its eigenbasis is of the form 



KWR F R W*K 



KWR y/foU v/p^J W V T 



1 

71 



-e-P 















(6.102) 
(6.103) 



e—+e-V/ 
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and the unitary matrix V equals rIw^K^TRWRqU, where T is the signature 
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1 











\ 







1 


















1 


















-1 


I 



(6.104) 



We have now all the necessary matrices to calculate the operator po — Fq in the Fq's eigenbasis. 
We obtain 



KWR (po-Fo)R W*K 



KWRqPqRIwW -KWRoy/pSUy/pdRlwitfT (6.105) 

KWPA 2 PW*K f -KWPAUAPW*K f T 

( Cosh(^)-Cos 2 (^) -iSin 2 (ji) \ 



e 2 



V 



2 j V2 





iSin 2 (fi) Sinh(%)-Sin 2 (%) j 



We are very closed to find the spectrum of po — Fq. We can denote by M the previous matrix. 
The eigenvalues of this matrix M are given by the roots of the polynomial P(x) = x 2 — Tr(M)jc + 
Det(M) which simply are 



x± = X - (tt(M) ± ^Tr(M) 2 -4Det(M)^j . 



(6.106) 



All this complicated construction was necessary to obtain the spectrum of the operator po — Fq. 
We victoriously end up with 

Spect(pQ-Fo) = i ^1 -e^ ± e ~ li yfl + -2e^Cos(^ . (6.107) 



This spectrum is not always positive (see Fig. 16.61) . 

Only in the regime of relatively large /i, the quantity 



e ^ y 1 + e ^ — 2e 2 Cos(n) is greater than 0. More precisely, 

Spect(p -F ) >0^/i >/!o~ 1.4386 

where /1q is the solution of the equation j ( 1 — — \j 1 + — 2e^~Cos(iA)^j = 



(6.108) 



In the regime jU > jlQ (see Fig. 16.71) . the optimal failure probability reaches the overall lower 
bound and we therefore get 



Q = F = Tr(Fo) = e~- 



(6.109) 
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The corresponding optimal measurement is moreover given by 



E = Z- l v ^(p -F ) v fp^- 1 (6.110) 

Ei = UE U 

E? = 1-Eq-UEqU. 



Note that for ji = jiQ, the POVM elements Eq and E\ have rank 1 since one eigenvalue of Pq — Fq 
vanishes. 
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We can remark here again that if we wanted to unambiguously discriminate 

Po = r O a X a l + \ioc){ia\) (6.111) 

and pi = - (| - a){-a| + | — ioc)(— ia\). (6.112) 
then we would find for jU > 0.7193 

Q = e-* i . (6.113) 



In the regime /i < /io where the operator po — Fq is not positive, we have to check the spectrum 
of the operator P^UP^. It is actually, as far as we know, not possible to calculate analytically its 
spectrum. Even if it is not really satisfying, we compute numerically the spectrum of P^UP^~. It 
turns out that it always has two eigenvalues of opposite sign in the regime /i < /lo- Consequently, 
we can write the operator P^UP^ in its eigenbasis {|0), |1)} as 

pj-UPf- = a|0)<0|-fc|l)(l|, fl,kR + . (6.114) 

And in virtue of Theorem 19, the optimal failure probability (see Fig. 16.81) for unambiguously 
discriminating the bit value mixed states is 

e opt = l-^ 7 (^(0|po|0)+a(l|p |l) + 2^|(0|po|l)|). (6.115) 
a + b 



Q 

1 

0.8 
0.6 
0.4 
0.2 



2 4 6 8 10' 

Figure 6.8: Optimal failure probability for USD of the bit value mixed states 

So far, no neat expression in terms of /i is known for this optimal failure probability <2 opt for 
M < Mo even if we do know its structure. This comes from the rather complicated form of the 
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states po and p\. As a final word, let us add that the optimal USD measurement is of form 



Eq = \x) (x\ 

Ei = UEqU with \x) 

E? = t-E - UE U 



/ e -/Arg({l|po|0)) \ 

Vl+a/b 
1 



V 








(6.116) 



/ 



even here also, we can note write them in term of the mean photon number /I. On the last graph 
16.91 we can show and compare the two optimal failure probabilities derived in this chapter. 




Figure 6.9: Comparison between the optimal failure probabilities for USD of the basis and the 
bit value mixed states 



This conclude the last chapter of this thesis. 



This last example might appear a bit unsatisfactory to the reader since no analytical expres- 
sion for P^UP^- is known. However this is exactly the contrary. During my work on Unambigu- 
ous State Discrimination, I was guided by the four density matrices presented in this chapter. 
They were my inspiration as well as my life ring. They are actually at the core of the derivation 
of the two classes of exact solutions and the numerous theorems derived in this thesis would not 
have been found without them. 
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Chapter 7 
Epilogue 



The main results of this thesis are, first, the two classes of exact solutions, second the reduction 
theorems, and finally the solution to unambiguous comparison of n pure states having some 
simple symmetry and the application of our results on USD to a BB84-type protocol. 

There are actually two directions for research in USD. The first path is of course the 
derivation of new solutions. The second is to find new applications of the already known 
solutions. In this thesis, we have tried to follow both paths. On one hand, we have derived new 
tools and new classes of exact solutions. On the other hand, we have given two examples of 
application for our tools. 

With respect to the newly developed tools, we have presented the notion of parallel addition 
p<yL~ l p\ m the context of unambiguous state discrimination. We have also shown the relevance 
of the two operators J ^/poPi ^/po and J ^fp[pQ^fp~\. We have finally provided two new classes 
of exact solutions as well as the three reduction theorems as we now discuss. 



The two classes of exact solutions derived in this thesis are the only two analytical solutions 
for unambiguous discrimination of two generic density matrices known so far. There now 
exist six analytical solutions for optimal unambiguous discrimination of quantum states. They 
correspond to the unambiguous discrimination of: 

1. Any set of linearly independent symmetric pure states |[T9l . 

2. Any pair of nonoverlapping mixed states^ such that the two operators po — oc J y/popiy/p~o 
and pi — li\/y/pipOy/p~i are positive semi-definite, and where a depends on the regime of the 

ratio ^ ^- [chapter 4]. Note that the case of 'Any pair of two pure states' solved by Jaeger and 



Any USD problem of two density matrices can be reduced to such a form according to Theorem 9. 
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Shimony lfT71 is included in this class of solutions. 

3. Any pair of geometrically uniform mixed states of rank two in a four-dimensional Hilbert 
space [chapter 5]. We find that only three options for the expression of the failure probability 
exist. First, if the operator po — y^^b is positive semi-definite, then the pair of density matrices 

falls in the first class of exact solutions. If this is not the case, either the operator P^UP^- has one 
positive and one negative eigenvalue or it has two eigenvalues of the same sign. In the former 
case, we can give the optimal failure probability in terms of the eigenvalues and eigenvectors of 
P^UP^. In the later case, no unambiguous discrimination is possible and the failure probability 
simply equals unity. 

4. A pure state and a density matrix with arbitrary a priori probabilities ll34l . 

5. Any pair of mixed states with one-dimensional kernel [|26l . 

6. Any pair of subspaces [|35l . 

Note that for the classes 2 and 3, we provide the optimal failure probability as well as the 
optimal measurement. Moreover, the solutions 4, 5 and 6 are reducible to some pure-state 
solutions. As we showed in this thesis, the reduction theorems and the solution for USD of two 
pure states are sufficient to derive those three solutions. 

The three reduction theorems allow us to reduce USD problems to simpler cases for which 
the solution might be known. This is the case, as we showed in chapter 3, for the unambiguous 
comparison of two pure states ll27l l28l l29l . the unambiguous comparison of n pure states 
having some simple symmetr]^, state filtering 031 l34l and the unambiguous discrimination 
of two subspaces lf35l . The reduction theorems also permit us to define a so-called standard 
USD problem. This problem is concerned with two density matrices of the same rank r in a 
2r-dimensional Hilbert space. This is proposed as a starting point for further investigations in 
unambiguous state discrimination in order to avoid trivial cases or unnecessary complexity. The 
reductions come from simple geometrical considerations and can be summarized as follows. 
With the first reduction theorem, we split off any common subspace between the supports of 
the two density matrices po and p\. Thanks to the second reduction theorem, we eliminate, if 
present, the part of the support of pi which is orthogonal to the support of po and vice versa. 
With the third reduction theorem, if two density matrices are block diagonal, we decompose 
the global USD problem into decoupled unambiguous discrimination tasks on each block. 
These three reduction theorems are also used to derive general theorems on unambiguous state 



2 n linearly independent pure states with equal a priori probabilities and equal and real overlaps. 
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discrimination. For example, the first reduction theorem is required to derive the two classes 
of exact solutions since the assumption of two density matrices without overlapping supports is 
made. 

With respect to the applications, we have used our new tools for the unambiguous compari- 
son of n pure states with a simple symmetry^ and to answer two crucial question^ related to the 
implementation of the Bennett-Brassard 1984 Quantum Key Distribution protocol. In fact we 
prove that the comparison of n linearly independent pure states with equal a priori probabilities 
and equal and real overlaps, a task related to the USD of two density matrices, can be reduced to 
n unambiguous discriminations of two pure states and can then be solved. The question to know 
whether any unambiguous comparison of pure states is always reducible to some pure state cases 
remains operJE With respect to the BB84-type protocol implemented with weak coherent pulses 
and a phase reference, we give the probability with which an eavesdropper can unambiguously 
distinguish the basis of the signal as well as the probability with which an eavesdropper can 
unambiguously determine which bit value is sent without being interested in the knowledge of 
the basis. 

Finally, as we discussed in chapter 5, a unified expression for the failure probability for the 
second class of exact solutions might be a pre-condition to find new solutions in unambiguous 
discrimination of two density matrices. Moreover new consequences of Theorem 18 should be 
investigated. 



3 « linearly independent pure states with equal a priori probabilities and equal and real overlaps. 

4 First 'With what probability can an eavesdropper unambiguously distinguish the basis of the signal?' and 
second 'With what probability can an eavesdropper unambiguously determine which bit value is sent without being 
interested in the knowledge of the basis?'. 

5 while the unambiguous comparison of mixed states is generally not reducible to some pure states case l28l 
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Chapter 8 
Appendix 



8.1 Appendix A 

Theorem 20 Theorem For any operator A, 

A*A\x) = 0^A t |x) = 0. (8.1) 



Proof We show this equivalence by proving separately the two implications. 
•<=] This direction is trivial. If A\x) = then A^Ajx) = 0. 

:] Here we make use of a fundamental theorem of linear algebra for any linear map A, the 
kernel of Af equals the orthogonal complement of the image of A that is to say Ker(A^) = 
7m (A) ^. Let us start with a vector \x) such that A'A|jc) = 0. A\x) is in the kernel of A^ so that 
A\x) is in Im(A)^. Moreover, by definition, A\x) is in Im(A). It implies that A\x) = 0. This 
completes the proof. ■ 



8.2 Appendix B 

Proof of Lemma 2 For any operator A, we can introduce a polar decomposition A = |A | V with 
|A| = VAA^ = V vAXAV^. Note that V is unitary and not necessarily unique, while vAA^ and 
VaU are unique and positive semi-definite. Moreover, since |A| might not have full rank, let us 
introduce the unitary transformation V' = ZV where Z is a unitary matrix of the form 




(8.2) 
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and T, a unitary matrix having support on ^J^i- From this remark, it follows that if A — \A\ V 
is a valid polar decomposition then A = \A\ V' is as well a valid polar decomposition. Indeed, 
A = \A\ V' = \A\ V and \A\ = V'Va^AV^ = V v / AUV t . 

We can now introduce a polar decomposition of A in the quantity Tr(AW) and find 

|Tr(AW)| = |Tr(|A|VW)| = |Tr(|A| 1/2 |A| 1/2 VW)|. (8.3) 

We denote X = \A\ [ ' 2 = and Y = \A\ l ^ 2 V W (W and V are both unitary matrices) and apply 
the Cauchy-Schwarz inequality (Theorem 2) to obtain 

|Tr(AW)| = iTr^rjl < ^Tr(|A|) ^Tr(WW\A\VW)) =Tr(|A|) . (8.4) 

Equality holds if and only if |A| 1/2 = fi\A\ l/2 VW, for some j3 £ C. This is possible if and only 
if fiVW = R, where R is of the same form than the unitary Z in Eqn. (I8.2I) . We can multiply each 
side with its adjoint and then find |/3 1 2 = 1. This implies that /3 = e~ l< ^ for some angle <j> so that 
we find the connection W = V^Re'^ . Since V comes from a polar decomposition of |A| and R 
is of the form of T, W* is a valid unitary for a polar decomposition of |A|. This completes the 
proof. ■ 



8.3 Appendix C 

Proof of Lemma 3 To complete the proof, we see two basic properties of the supports of two 
positive semi-definite matrices M and iV 

-Ymn C rY M , (8.5) 



y M c y M+N . (8.6) 

The first ingredient is to see that A : B is Hermitian. Indeed, we can write 

A(A + B)- l B = A(A + B)-\B + A-A) (8.7) 

= A(A + By 1 (A + B)-A(A + By 1 A. (8.8) 

Let us underline that A (A + B)~ l (A + B) = AYl r y A+B = A since ,y A C S^a+b- Similarly (A + 
B){A + BY l A = Uy A+B A=A. As a result, 

A(A + BY l B = A-A{A + By l A (8.9) 

= (A + B)(A + By 1 A-A(A + By 1 A (8.10) 

= A(A + B)- l B. (8.11) 
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Now we can prove that S^a-.b C S^a H =5^. Indeed ^a(a+b)-is C and ^b(a+s)-ia c ^s- 
Since A (A + 5) = B(A + B)~ l A, it follows that ^ A:S C J^H^. 

The last step is to prove that 5? \ PI <5*b C S^a.b- To do so, let x be in ,5^ fl 5?b and find a vector 
y E S?a U such that (A : 2?)y = x. Actually, such a y is given by (A -1 + B 1 )x. Indeed 

(A:B)y = A(A + B)~ l B(A~ 1 +B~ 1 )x (8.12) 
= B(A + By l AA~ l +A(A + By 1 BB~ 1 x (8.13) 
= B(A + B)- l x + A(A + B)- l x (8.14) 



since, Wx E S^a H AA *x = jc and BB l x = x. Finally we can write (A : B)y = (B+A)(A + 
B)~ l x = x since x E S^a H ^b C =^+e- These completes the proof. ■ 
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